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The  purpose  of  this  effort  was  to  demonstrate  that 
moving-bank  multiple  model  adaptive  estimation  algorithms 
could  be  applied  to  a  realistic  and  practical  control 
problem.  Moving-bank  multiple  model  adaptive  estimation/ 
control  is  an  atempt  to  create  from  full-bank  multiple  model 
adptive  estimation/control  an  adaptive  e s t ima t  ion/ con t ro 1 
technique  which  maintains  its  desirable  on-line  parameter 
adaptation,  but  which  reduces  the  intensive  computational 
loading  which  makes  the  full-bank  method  impractical.  The 
results  of  this  thesis  do  indicate  that  the  moving-bank 
multiple  model  adaptive  algorithm  can  be  successful  in  this 
goal. 
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for  the  help  he  gave  me  at  every  step  of  this  effort,  and 
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would  also  like  to  thank  my  classmates  for  thier 
encouragement,  especially  S.  Eckert  for  helping  me  procure 
valuable  computer  resources. 
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Abstract 

This  investigation  applies  moving-bank  multiple  model 
adaptive  estimation/control  algorithms  to  the  control  of  a 
realistic  model  of  a  large  flexible  spacecraft.  Moving-bank 
multiple  model  adaptive  estimation  differs  from  conventional 
(full-bank)  multiple  model  adaptive  estimtion  in  that  a 
substantially  reduced  number  of  elemental  filters  is 
required  for  the  moving-bank  estimator  (9  vs.  100  for  the 
system  modeled  in  this  thesis).  The  positions  in  parameter 
space  that  the  reduced  number  of  elemental  filters  occupy 
are  dynamically  re-declared;  i.e.,  the  moving-bank  slides 
about  the  parameter  space  in  search  of  the  true  parameter 
point. 

^Critical  to  the  performance  of  the  moving-bank  multiple 
model  adaptive  estimator  is  the  decision  logic  used  to 
determine  which  elemental  filters  are  implemented  in  the 
bank,  and  when  to  change  this  decision.  The  decision  logics 
discussed  focus  on  three  situations:  initial  acquisition  of 
the  unknown  parameter  values,  through  reducing  bank 
discretization;  tracking  the  unknown  parameter  values, 
through  bank  movement;  and  r e a c qu i s i t i on  of  the  unknown 
parameters  following  a  large  jump  change  in  their  values, 
through  expanding  bank  discretization.  Ambiguity  function 
analysis  is  used  to  predict  performance  in  these  situations. 

v  i  i 


-  The  system  to  be  controlled  is  a  simplified  model  of  a 
large  scale  space  structure.  Its  equations  of  motion  are 
developed  and  placed  in  state  space  form,  the  states  being 
the  positions  and  velocities  of  the  rigid  body  mode  and  the 
second  and  fourth  bending  modes.  The  state  space  matrices 
describing  the  system  are  computed  based  on  nominal  values 
for  all  physical  parameters  with  the  exception  of  the  mass 
density  of  the  structure  arms  and  their  modulus  of 
elasticity.  These  two  parameters  are  allowed  to  vary  in 
discrete  steps,  estaplishing  the  parameter  space.  It  is 
then  attempted  to  control  the  states  to  the  quiesent  state, 
using  moving-bank  multiple  model  adaptive  algorithms. 

The  results  indicate  that,  although  the  system  under 
study  did  not  have  a  great  need  for  adaptive  estimation  and 
control,  the  multiple  model  adaptive  estimator  performs 
essentially  identically  to  a  single  filter  artificially 
knowledgeable  of  the  uncertain  parameter  values.  In 
addition  changing  bank  discretization  for  the  initial 
parameter  aquisition  phase  speeded  acquisition.  However,  th 
bank  was  unable  to  expand  following  a  jump  change  in  the 
uncertain  parameter  values,  in  order  to  restart  the 
acquisition  phasep  the  bank  tracked  the  jump  change  through 
movement  alone.  Ambiguity  function  analysis  proved  to  be  an 
excellent  predictor  of  bank  performance,  and  should  be  used 
as  a  design  tool. 


v  i  i  i 


MOVING-BANK  MULTIPLE  MODEL  ADAPTIVE  ALGORITHMS 
APPLIED  TO  FLEXIBLE  SPACECRAFT  CONTROL 


I .  Introdac  t  ion 

A  common  problem  in  estimation  and  control  problems  is 
the  uncertainty  of  parameters  used  in  the  design  of  the 
system  model  and  embedded  in  the  estimator  or  controller. 
These  parameters  can  be  initially  unknown  but  unchanging, 
or  they  can  vary  slowly,  or  they  can  undergo  abrupt  changes 
as  in  the  case  of  partial  system  failure.  If  these 
parameters  vary  enough,  it  may  be  necessary  to  estimate 
their  values  along  with  the  system  state,  and  adapt  the 
estimator,  controller,  or  both,  to  incorporate  the  current 
value  of  the  uncertain  parameters.  This  requirement  is  not 
based  soley  upon  the  magnitude  of  the  parameter  variation 
but  also  upon  the  sensitivity  of  the  system  to  the 
variation.  Some  parameters  may  vary  widely  with  no 
degradation  in  system  performance,  while  a  small  variation 
in  another  parameter  may  cause  the  system  to  become 
unstable.  This  thesis  expands  the  exploration  of  one  method 
for  accomplishing  this,  known  as  moving-bank  multiple  model 
adaptive  estimation. 

I . 1  Background 

In  a  large  class  of  problems,  which  can  be  modeled  as 
linear  stochastic  systems  with  uncertain  parameters 
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The  convergence  of  the  MMAE  algorithm  estimate  to  the 
true  state  has  been  shown  for  the  case  of  uncertain  but 
constant  parameters  [6,  7].  Convergence  has  not  been 
shown  for  varying  parameters;  however,  promising  results 
have  been  obtained  for  an  ad  hoc  approach  where  constant 
parameters  are  assumed  for  algorithm  design,  but  the 
computed  conditional  probabilities  are  lower  bounded  to 
prevent  'locking  onto'  a  single  parameter  value  [5,8,9]. 

If  the  possible  parameter  time  variations  are  modeled, 
the  optimal  state  estimate  can  be  obtained  as  the  weighted 
sum  of  the  estimates  produced  by  filters  matched  to  all 
possible  parameter  time  histories  [5].  This  approach 
is,  however,  impractical  as  K1  elemental  filters  at  sample 
time  t^  would  be  required.  Even  when  the  parameter  temporal 
variation  can  be  well  modeled  as  a  Harkov  process,  it  has 
been  shown  that  the  number  of  elemental  estimators  required 
would  be  K  [5,10].  Thus  if  only  2  uncertain  parameters 
were  each  discretized  to  10  values,  K  would  be  equal  to  100, 
and  K  =  10,000:  still  impractical  for  implementation  [5]. 

For  control  applications,  the  state  estimate  obtained 
from  a  multiple  model  adaptive  estimator  can  be 
premultiplied  by  a  controller  gain  established  via  forced 
equivalence  design  [12].  The  gain  may  be  based  on  a 
single  nominal  parameter  value,  or  can  be  evaluated  using 
the  estimated  values  of  the  uncertain  parameters  provided 
by  the  estimator.  If  a  separate  controller  gain  is 
associated  with  each  elemental  filter  in  the  bank,  the 


control  can  be  produced  as  the  probabilistically  weighted 
average  of  the  elemental  f i 1 t e r / cont r o 1 1 er  outputs,  in  the 
same  way  the  state  estimate  is  obtained  in  a  multiple  model 
adaptive  estimator.  This  is  known  as  multiple  model 
adaptive  control  (MMAC)  [12:253,11].  An  alternative,  using 
MAP  criteria,  is  to  select  that  f i 1 1 e r / c on t r o 1 1  e r  output 
associated  with  the  highest  weighting  probability  [11]. 

MMAE  has  successfully  been  applied  to  several  practical 
problems.  The  tracking  of  maneuvering  targets  has  been 
shown  to  lend  itself  to  this  approach  [13,14,15,16,17]. 

Other  applications  have  been  demonstrated  in  flight  control 
[11],  multiple  hypotheses  testing  [18],  detection  of 
incidents  on  freeways  [19],  adaptive  deconvolution  of 
seismic  signals  [20],  and  problems  in  which  initial 
uncertainties  are  so  large  that  nonadaptive  extended  Kalman 
filters  diverge  [21,22]. 

In  implementations  where  there  are  a  large  number  of 
uncertain  parameters,  or  where  the  uncertain  parameters  can 
take  on  a  large  number  of  discrete  values,  the  storage 
required  for  the  elemental  filters,  and  the  computational 
loading  required  become  unwieldy.  Consider,  if  there  are 
only  2  uncertain  parameters,  each  of  which  can  assume  10 
discrete  values,  the  number  of  elemental  filters  is 
10  -  100;  this  may  be  i m p 1 e m e n t a b  1  e  but  it  can  be  seen 

that  the  number  of  filters  required  will  grow  quickly  with 
additional  uncertain  parameters  or  finer  levels  of 
discretization.  Approaches  to  reducing  the  computational 
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burden  include:  use  of  Markov  models  for  parameter 
variation  [5,10,23];  'pruning'  and  'merging*  of  branches  in 
a  'tree'  of  possible  parameter  time  histories  [24,25], 
hierarchical  structuring  [26],  and  dynamic  co ar s e- to- f ine r 
rediscretization  [27]. 

A  method  proposed  by  Hentz  and  Maybeck  [28,29] 
involves  implementing  only  a  smaller  number  of  elemental 
filters  selected  from  the  total  bank.  This  'moving-bank'  is 
adjusted  by  replacing  filters  currently  implemented  with 
others  from  the  larger  bank  in  response  to  the  changes  in 
the  weighting  probabilities  and  filter  residuals.  In  the 
above  example,  only  the  filters  with  the  three  closest 
discrete  values  to  each  of  the  estimated  parameter  values 
might  be  implemented.  This  would  result  in  only  nine 
filters  being  implemented,  and  as  the  estimate  of  the 
parameters  changed,  the  implemented  bank  of  filters  would 
'move';  see  Figure  1-1.  This  approach  was  shown  to  achieve 
performance  essentially  equivalent  to  that  of  a  single 
Kalman  filter  artificially  knowledgeable  of  true  parameter 
values,  for  a  simple  but  physically  motivated  example  of  a 
single  input/single  output  second  order  system  with 
uncertain  damping  ratio,  and  undamped  natural  frequency. 

This  corresponds,  for  example,  to  a  bending  mode  in  an 
aerospace  vehicle  [28:15,57-101;  29:17-27].  The  algorithms 
used  for  implementing  the  moving-bank  are  developed  in 
detail  in  Chapter  II. 


□  elemental 
Kalman 
f i I  ter 

X  true 

uncertain 
parameter 
val  ue 


Figure  1-1.  Moving-Bank  Multiple  Model  Adaptive  Estimator. 


I • 2  Prob  lea 

Full-scale  MMAE  requires  too  large  a  computational 
load  to  be  practical  for  most  applications  requiring 
adaptive  e s t im a t ion/ c on t r o 1.  Moving-bank  MMAE  has  been 
shown  to  be  feasible  for  a  simple  but  physically  motivated 
application.  However,  it  still  needs  to  be  demonstrated 
that  moving-bank  MMAE  can  be  useful  in  a  more  complicated 
and  realistic  application  requiring  adaptive  estimation/ 


control . 


The  purpose  of  this  effort  is  to  apply  moving-bank  MMAE 
to  a  practical  application  with  a  higher  level  of  complexity 
and  farther  evaluate  its  performance  potential.  In  addition 
alternative  decision  logics  for  chosing  which  of  the 
elemental  filters  in  the  parameter  space  to  implement  in  the 
moving-bank  will  be  evaluated. 

1 . 3  Scone 

A  moving-bank  multiple  model  adaptive  estimator/ 
controller  is  evaluated  for  a  realistic  application.  A 
simplified  model  of  a  large  flexible  spacecraft,  consisting 
of  a  central  rigid  hup  with  four  radiating  flexible  arms,  is 
used.  The  flexible  arms  are  representative  of  the  flexible 
appendages  such  as  solar  panels  or  antennas  attached  to 
actual  spacecraft.  The  spacecraft  is  described  in  terms  of 
its  physical  parameters  (mass,  height,  length  of  the  arms, 
etc.),  with  two  uncertain  parameters  (mass  density  of  the 
arms,  and  the  modulus  of  elasticity  of  the  arms).  The 
uncertain  parameters  are  discretized  into  10  values  each,  to 
provide  a  ten  by  ten  (100  point)  parameter  space.  Three 
point  actuators  (representative  of  the  pulse  rocket  motors 
used  on  actual  spacecraft)  provide  the  control  input,  and 
five  position  sensors  along  with  five  velocity  sensors 
provide  the  measurements  for  updating  the  state  estimates. 
The  dynamics  and  measurement  noise  characteristics  are 
assumed  known  with  no  uncertainty,  and  are  modeled  as  white 
Gaussian  processes  with  strengths  Q  and  R  respectively. 
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I  .  4  £££££££& 

The  perfomance  of  the  moving-bank  mnltiple  model 
adaptive  eitimator/controller,  its  ability  to  estimate  both 
the  states  and  parameters  as  well  as  to  apply  adequate 
control  to  the  true  system  is  compared  to  a  benchmark  of  a 
single  Kalman  filter/controller  based  on  (artificial) 
perfect  knowledge  of  the  trne  parameter  value.  Estimator 
performance  is  evaluated  under  the  following  conditions: 

a.  The  true  parameter  value  is  constant  and  exactly 
equal  to  one  of  the  discretized  parameter  values  within  the 
initial  conditions  chosen  for  the  bank. 

b.  The  true  parameter  value  is  constant  and  lies 
outside  the  initial  conditions  chosen  for  the  bank  (but 
within  the  overall  parameter  space). 

c.  The  true  parameter  value  undergoes  jumps  in  value 
from  one  of  the  discretized  parameter  values  to  another, 
the  values  chosen  so  that  the  new  true  parameter  point  is 
outside  the  area  covered  by  the  moving  bank. 

True  parameter  values  which  are  not  exactly  equal  to  one  of 
the  discrete  parameter  values  must  be  investigated  before  a 
full  evaluation  of  moving-bank  multiple  model  adaptive 
estimation/control  can  be  done  [28];  however,  due  to  time 
constraints,  they  are  not  treated  in  this  thesis. 

The  effect  on  performance  and  computational  burden  of 
various  algorithms  that  can  be  used  for  initial  acquisition 
of  the  true  parameter  values,  and  for  identification  of  when 


a  jump  change  has  taken  place,  are  of  interest.  Specific 
algorithms  developed  in  Chapter  II  vill  be  tested  to 
determine  their  effect  on  estimator  performance. 

Comparison  of  contoller  design  strategies,  contrasting 
a  single  fixed— gain  controller,  a  single  change ab 1 e— ga in 
controller,  and  a  moving-bank  multiple  model  adaptive 
controller,  vill  also  be  accomplished. 


I  •  5  Or g.AftiiAiiafi 

The  remaining  sections  of  this  thesis  are  organized  as 
follows.  Chapter  II  develops  the  algorithms  nsed  for  both 
full-scale  MMAE  and  moving-bank  MMAE.  Chapter  III  presents 
the  flexible  spacecraft  structure  model.  Chapter  IV  details 
the  siaulat' m  used  for  the  moving-bank  e s t i m a to r / cont r o 1 1 e r 
performance  evaluation.  Chapter  V  contains  an  analysis  of 
the  results  of  the  evaluations.  And  Chapter  VI  presents 
conclusions  and  recommendations. 


II.  ALGORITHM  DEVELOPMENT 


1 1 .  i .  IaUfl4ji£tian 

This  chapter  details  the  algorithms  for  two  of  the 
adaptive  estimation  techniques  introduced  in  the  first 
chapter.  First  full-scale  Bayesian  Multiple  Model  Adaptive 
Estimation  is  developed,  then  the  modifications  necessary 
for  the  moving-bank  case  are  discussed. 

1 1.2.  Bave s ian  Est imat ion  Algor ithm  Deve lonment 

Full-scale  Bayesian  Multiple  Model  Adaptive  Estimation 
algorithms  are  presented  in  this  section.  For  a  more 
complete  development  of  concepts,  see  reference  [5:129-136]. 

Let  the  system  under  consideration  be  described  by  the 
following : 

A(ti  +  1)  =  ♦<ti  +  1.ti)i<ti)  +  Bd(ti)a(ti)  +  GjUpw^ti) 

iUj)  =  +  x.(  t  i ) 

IV 

where,  letting  the  denote  random  process: 


x  ( t  £ )  : 

state 

vector 

t  i  +  i ,  t  i  )  : 

state 

transiton  matrix 

ft  ( t  j )  : 

known 

input  vector 

®d ( t  i )  : 

control  input  matrix 

*d(ti>: 

wh  i  t  e 

Gaussian  dynamics  noise  vector 

noise 

input  matrix 
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measurement  vector 


i 


i( t4)  : 

H(ti)  : 

i<tA)  : 

a# 


measurement  matrix 

white  Gaussian  measurement  noise  vector 


and  the  following  statistics  apply: 

Etwd(ti)}  =  0 

E{v(ti)}  =  0 
E{£(ti)vT(tj)}  = 

where  5  —  is  the  Kronecker  delta  function.  It  is  also 
assumed  that  s.(tg),  w<j(ti),  and  v(t^)  are  independent  for 
all  tj. 

Now,  let  jj.  be  the  uncertain  p-dimens  ional  parameter 

a* 

vector  which  is  an  element  of  A,  where  A  is  a  subset  of  Rp. 

£  may  be  uncertain  but  constant,  it  may  be  slowly  varying, 
or  it  may  under  go  jump  changes.  The  parameter  vector  a  can 
affect  any  or  all  of  the  following:  6,  B  d ,  G  d ,  Qd,  B,  and  R. 
The  Bayesian  estimator  computes  the  following  conditional 
probability  density  function: 

£i< t .ii|( tt> ■  fi< tt) I tj) 

where  Z(t^)  is  a  vector  of  measurements  from  tg  to  t^, 

*0 

Z(t£)  =  [*T(ti)  .^(t-.j) _ ,zT(t0)]T 
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The  second  tern  on  the  right  side  of  Equation  (II-l) 


can  be  further  evaluated: 


fg.|Z<  tA) 

«N»  <*■*  1 


fA.i(ti)  Igttj.i)  (A*ii^i-l) 
fi(ti)lg(ti_1)(Ai,^i-l) 

fi(  ti)  U.g(  t  i_x  >  (Ai  I  A'Zi-i>  filg(  t  i.j)  (*lli-l) 
ifz(t.)la,Z(ti_,)(£ilA'£i-l)fa|Z(ti_1)<A|£i-l)iA 

Am  1  M  IX  M  M  AX 


(II-2) 


Conceptually,  Equation  (II-2)  can  now  be  solved 
recursively,  starting  from  an  a  priori  probability  density 

of  f&(a),  since  f  *(  t  .  )  I  a,Z(  t  {i < S-i  I  A*^i-1 )  is  G»BSsian  with  a 
mean  of  H(ti)x(ti”)  and  covariance  [ H( t P( t j-) HT( t A ) +  R( t ]  , 
where  j(tj”)  and  PCtj-)  >re  the  conditional  mean  and 
covariance  respectively  of  x(tj)  just  prior  to  the 
measurement  at  t-,  assuming  a  particular  realization  a.  of  j». 

Using  the  conditional  mean,  the  estimate  of  i.(t^)  can 
be  gene  rated  as: 

EtiUilIZUi)  =  ZiJ  =  I  ifI(t.)|z(ti)(i|Zi>4i 

M#  —09  »  1  -  * 

=  I-a,x[/Afx(ti),a|Z(ti)(A*Al2i)daUi  (II-3) 
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(  d  I  Z  )  cal  dj. 


_  f  f*  x-  *  x  Cxi  a.Z  .  )  cxl  *  .  ■  2{t 

-  Jav  —  x(t;)ici.Z(t)  -  — *  al  ^  ^ T  i 


( 5. 1  Z  ;  >  c  a 


(II  4> 


where  the  term  in  brackets  is  the  estimate  of  s.(t^)  based  on  a 


particular  value  of  the  parameter  vector.  This  would  be  the 


output  of  the  Kalman  filter  based  on  a  particular  parameter 


value.  When  a.  is  continuous  over  A,  this  would  require  an 


infinite  number  of  filters  in  the  bank.  In  the  continuously 


distributed  parameter  case,  what  is  typically  done  is  to 


discretize  the  parameter  space,  yielding  a  finite  number  of 


filters.  The  integrals  over  A  in  Equations  (II-2)  through 


(II-4)  then  become  summations:  letting  p^(t^)  be  defined 
as  the  probability  of  the  k^  elemental  filter  being 


correct,  conditioned  on  the  measurement  history,  (II-2)  and 


(II-4)  are  replaced  by: 


fz(ti)  I  a,Z(ti_1)  Ak'^i-l)Pk(ti-l) 


Pk(ti)  = 


(II-5) 


^  f  A  ( t  i )  la.,Z(ti_1)  ^  —  i  ^  —  j  '  — i-l^Pj  ^  ti-l' 

j=l  -  ~  ~ 


x  <  t  i+ )  =  EtxUi)  I  Z  ( t  i  >  =  Zt} 


=  \  Ak(ti  +  )pk(ti) 


( I 1-6 ) 


where  a  [  Aj » &2  *  •••  *  -K  ^  an<*  £k(ti+)  is  the  mean 


onditioned  on  a  =  ak  and  Z(t^)  =  Z^,  i.e.  the  output  of  the 


kc  Kalman  filter  in  the  bank,  based  on  the  assumption  that 


ji  =  ak.  Pictorially  the  algorithm  appears  as  in 


F igure  II-l. 
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Figure  II-l.  Multiple  Model  Filtering  Algorithm 


The  probability  weighting  factors  for  each  Kalman 
filter  are  calculated  from  Equation  (II-5),  where: 


fz(ti)  la.Kti.!)  (^i**k'^i-i: 


(2n)m/2|Ak( t±) I (1/7> 


tp[— (1/2)  rkT  (ti)Aj,-1(ti)  r.k  ( t  £  )  ]  (II-7) 


Ak(ti}  =  Hk(ti)Pk(ti')Hki(ti)  +  Rj.Ctj) 


rk(ti>  =  z.£-  Hk  (  t  i  )  xk  ( t  i  ) 


m  =  dimension  of  z  (the  number  of  measurements) 


Both  the  residual  covariance  Ak(t^)  and  the  residual  rk(tj) 
itself  are  readily  available  from  the  k**1  elemental  filter. 


The  estimate  of  the  parameter  and  the  covariance  of  the 
parameter  are  given  by: 


-  f 


aUi)  =  E  { a  I  Z  ( t  i )  =  Zi)  =  J  a  f g  ( t  .  )  (  a  I  Z  £ )  ia 

'v  —•  qo  33*  *5  1 

K 


=  I  a  [  ^  Pjc(ti)6(a-  ak 


)  1  4a 


-00  k  =  1 

K 

^  akpk(ti) 


k=l 


and 


E { [ a  -  a ( t  £ ) ] [a  -  £( t i ) ] T I Z ( t  i )  =  Z = 


^  [ak  -  attj^Jlta^-  aC^ 
k=l 


)]  Pk<*i> 


The  covariance  of  the  state  estimate  is  given  by: 
PU^)  =  EUxU^  -  £(ti+)l  [xttj)  -  £(ti+)]T|Z(ti)  -  ; 

W 

=  J  [x  -  x(t.  +  )Hi  -  i(ti  +  )]Tfi(t  i|2(t  )  (xlz 
—  00  *>' 

K 

=  ^  Pk(ti>(  [ x  -  x(ti+)][x  -  x(ti+)]T 
k=  1 

f  ^.(t  lla,|(t 


=  ^  Pfc  (  t  i )  {Pk  ( t  i+  )  +  [x  ( t  k  + )  -  x(ti+)] 

Cxk( t i  +  )  -  x( t  A  +  ) ]  T) 


k=  1 


where  Pk(t^+)  *s  covariance  of  the  state  estimate 


th 


( I I- 8 ) 


( 1 1-9 ) 


i)dx 


i )  dx } 

(11-10) 
of  the 


elemental  filter 


It  is  important  to  note  that  it  is  not  necessary  to 


compute  the  covariance  of  the  parameter  estimate  and  the 
covariance  of  the  state  estimate  online.  They  are  used  for 
the  evaluation  of  estimator  performance  daring  estimator 
design  and  tuning. 

1 1. 2.1.  Elemental  liliaaior.  The  elemental 

estimators  are  each  constant  gain- Kalman  filters,  the  gains 
of  which  are  dependent  upon  the  associated  point  in  the 
parameter  space.  Thus,  each  is  designed  on  the  basis  of  a 
time  invariant  system  model,  and  a  stationary  noise  model, 
and  the  initial  transient  in  the  resulting  filter  is 
ignored . 

The  k1*1  elemental  filter  estimate  of  x(t),  denoted  by 
xk(t),  is  propagated  from  just  after  the  last  measurement  to 
just  before  the  next  measurement  by: 

xk(ti_)  =  *kxk(ti_1+)  +  Bdku(ti_1)  (11-11) 

and  is  updated  by: 

2k(t  i  +  )  =  xk(ti")  +  Kk  [  ^  (  t  ±  )  -  Hkxk(ti")l  (11-12) 

In  these  equations,  the  subscript  'k'  indentifies  each 
vector  and  matrix  as  being  specifically  based  on  the  k**1 
point  in  the  parameter  space. 


II. 2. 2.  Convergence.  The  adaptive  estimation 
algorithm  just  developed  is  optimal  under  the  conditions 
that  the  discretized  parameter  space  is  a  true  model  of  the 
physical  parameters,  and  the  true  parameter  vector  remains 
at  a  constant,  but  unknown,  value.  Under  these  conditions 
the  algorithm  will  converge  to  the  true  parameter  value 
[ 6  ]  ,  i. e . , 

lim  p^(tj)  =  0  for  a  ^  a^ 

1  -*  CD  ~ 

lim  p^(t^)  =  1  for  a.  =  a^ 

i  -*  CD  — 

When  the  true  parameter  space  is  in  fact  continuous  and  the 
true  parameters  lie  somewhere  between  the  discretized 
points,  the  algorithm  converges  to  the  single  discrete 
parameter  point  that  is  'nearest',  as  defined  in  [6],  to  the 
true  parameters. 

These  results  were  extended  by  Dasgupta  and  Westphal 
[30]  for  the  case  of  unknown  biases  in  the  measurement 
process,  (E(v(t^)}  =  mv(t^>,  where  my(t^)  can  be  affected  by 

the  parameter  vector).  Under  these  conditions  the  algorithm 
may  converge  to  a  parameter  point  that  is  not  close  to  the 
true  value,  and  erroneous  estimates  may  result. 

It  is  important  to  remember  that  the  algorithm  is  able 
to  identify  the  closest  parameter  point  only  by  observing 
which  elemental  filter  consistently  has  the  smallest  value 
of  Xfc  the  smallest  residual  relative  to  the  k 

filter-predicted  residual  squared  value.  If  the  residuals 
are  of  the  same  size.  Equation  (II-7)  shows  that  the  filter 
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with  the  smallest  value  of  |A^|  will  be  identified  as  the 
correct  filter.  Since  I I  is  independent  of  both  the 
residuals  and  the  elemental  filter's  correctness,  if 
pseudonoise  is  added  during  the  filter  tuning  process,  to 
account  for  model  inadequacies,  it  must  not  be  so  strong  as 
to  mask  the  cor re c tne s s / inc o r r e c tne s s  of  the  elemental 
filters.  Such  strong  pseudonoise  may  well  allow  the 
adaptive  filter  to  converge  to  an  erroneous  parameter 
value  [30] . 

As  noted  in  the  introduction,  no  satisfactory 
theoretical  convergence  results  are  available  for  more 
general  conditions,  such  as  slowly  varying  parameters, 
although  empirical  information  suggests  convergence.  The 
most  sucessful  approach  used  to  prevent  the  algorithm  from 
locking  onto  one  elemental  filter  before  the  parameters  have 
varied  significantly  has  been  to  lower  bound  the  p^'s  to 
prevent  them  from  converging  to  zero  [5,10]. 

II. 2. 3.  Contro 1.  There  are  several  'assumed 

certainty  equivalence  design'  [12:241]  approaches  to 
controlling  systems  with  uncertain  parameters  that  can  be 
used  with  a  multiple  model  adaptive  estimator.  In  the 
first,  the  estimator  provides  only  a  state  vector  estimate 
to  a  fixed-gain  controller  robustified  around  a  nominal 
value  of  the  uncertain  parameter  vector,  inom.  The 
controller  design  method  and  controller  gains  are 
independent  of  the  adaptive  nature  of  the  estimator. 


The  controller  algorithm  is  of  the  form: 

a<ti>  =  -Gc*[ti'Ano«li(ti+)  (11-13) 

or  the  steady-state  constant-gain  version: 

a<ti>  =  'Gc*Uno-li(ti+>  (II-14) 

A  second  approach  is  to  have  the  estimator  provide  both 
an  estimate  of  the  state  vector  and  an  estimate  of  the 
uncertain  parameter  vector  (Equation  II-8)  to  the 
controller.  The  controller  gain  then  becomes  dependent  on 
the  parameter  vector  estimate: 

a(  t  i)  *  -Gc*[ti,a(ti‘)]x(ti  +  )  (11-15) 

or  the  steady-state  c ons t ant- g a  in  version: 

a( t  A)  *  -Gc*t&(ti“)]£(ti  +  )  (11-16) 

where  a.(t|~)  (as  generated  at  the  previous  sample  time)  is 
used  instead  of  a(t^+)  to  reduce  computational  delay.  It  is 
important  that  the  control  input  be  applied  as  close  to  the 
start  of  the  sample  period  as  possible,  as  the  filter 
propagation  equation  is  written  on  the  assumption  that  the 
control  input  is  present  for  the  entire  sample  period. 

A  third  approach,  referred  to  as  multiple  model 
adaptive  control,  [12:253],  is  to  form  K  elemental 
controllers,  each  associated  with  one  of  the  elemental 
estimators  in  the  parameter  space.  Then  the  final  control 
input  becomes  the  probabilistically  weighted  average  of  the 


individual  controller  results  in  the  sane  manner  as  the 


state  estimate  was  achieved: 

*k(ti>  =  -Gc*[ti^k]ik(ti  +  )  (11-17) 

K 

»<t4)  =  1  P(ti)ak(ti)  (11-18) 

k=i 

For  this  thesis,  the  controllers  to  be  used  will  be 
linear,  quadratic  cost.  (LQ)  full-state  feedback  optimal 
deterministic  controllers,  designed  to  regulate  the  system 
to  the  quiescent  state,  incascade  with  the  state  estimator 
(invoking  assumed  certainty  equivalence).  The  basic 
structure  of  all  the  controllers  described  above  will  be 
similar,  where  the  gain  matrices  will  be  dependent  upon  a 


particular  parameter  vector  value. 

If  we  are  given  the  stochastic  system  [28:33-35]: 


x(t)  =  Fj,(t)  +  Bs  +  Gw(t) 


(11-19) 


where 


*,V 

'»> 

W! 


E  ( w  ( t )  )  =  0 


E { w ( t ) w( t+r ) }  =  Q8(x) 


and  the  quadratic  cost  function  to  be  minimized  is: 


J  =  E{  J 


(1/2) (£T( t)Wxx(t)  +  ST< t)Wuq( t) ]dt}  (11-20) 
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B(t. tA)  =  /  *(t. 

J  t{ 


r )  B  dT 


3d  "  B(ti+l*ti) 


♦(t2»ti)  is  the  state  transition  matrix  from  t2  to  t2  and: 

♦  =  ♦<ti+1,ti) 


Note  that  Equation  (11-21)  is  also  the  solution  to  the 
deterministic  LQ  optimal  control  problem  with  no  driving 
noise  g(t),  and  that  if  the  assumption  of  fall  state  access 
is  replaced  by  noise-corrupted  measurements  being  available, 
then  g(t^)  in  Equation  (11-21)  is  replaced  by  the  state 
estimates  i(t  generated  by  a  Kalman  Filter.  This  assumed 
equivalence  is  valid  if  all  system  parameters  are  known 
perfectly.  Therefore  for  the  uncertain  parameter  case  Gc* 
should  be  a  function  of  the  uncertain  parameter  vector  .a, 
and  'forced*  or  'assumed'  certainty  equivalence  design  is 
used  for  synthesis  [12:241]. 


1 1 . 3  .  Mfivi.ai.~Ba.filt  Al.gfiri.ilim  Deve  lopment 

As  the  number  of  uncertain  parameters  to  be  estimated 
grows  larger  and  the  discretization  of  the  parameter  space 
becomes  finer,  the  MMAE  algorithm  becomes  computationally 
impractical  for  real  time  applications.  Maybeck  and  Hentz 
[28,29]  have  demonstrated  that  it  is  feasible  to  begin 


with  the  entire  bank  of  Salman  filters  used  for  MMAE  but 


only  compute  the  state  estimate  and  parameter  estimate  for 
those  filters  'closest*  to  the  current  estimate  of  the 
parameter  vector.  The  p^’s  of  the  unimplemented  filter  are 
inherently  set  to  0,  and  all  of  the  probability  weighting  is 
distributed  amongst  the  implemented  filters.  As  the 
parameter  vector  estimate  changes,  the  implemented  bank 
'slides*  within  the  larger  bank:  those  implemented  filters 
'farthest*  from  the  current  estimate  of  the  parameter  vector 
are  dropped  from  the  bank  (no  longer  computed),  and  new 
filters  'closer*  to  the  current  parameter  vector  estimate 
are  implemented  instead.  It  is  also  possible  for  the  moving 
bank  to  change  discretization  level,  i.e.,  the  filters 
implemented  need  not  be  adjacent  in  the  full  bank.  Naybeck 
and  Hentz  examined  changing  discretization  level,  starting 
with  a  coarse  discretization  during  an  acquisition 
period,  then  changing  to  a  fine  discretization  once 
parameter  acquisition  was  achieved  [28,29]. 

II. 3.1.  ls.AAJl.iS.  4  A v e_XAAS.-  The  outputs  of  each 

elemental  filter,  £j(t^+),  are  put  through  a  weighted 
average  in  the  same  manner  as  in  Equation  (II-6),  except 
that  instead  of  summing  over  the  full  set  of  K  filters,  only 
those  implemented  in  the  moving  bank  are  summed.  Thus  if  J 


filters  are  implemented.  Equation  ( 1 1  —  6  )  becomes: 


£(ti+) 


=  \ ^ijtti+lp  j  <t  £) 

j  =  1 


(11-24) 


Similarly.  Equation  ( 1 1  —  5 )  describing  the  Pj(t^)’s  becomes: 


Pj(t.)  = 


f  j  (i(  tt)  )pj  ( ti_1! 


(11-25) 


))Pk(ti-1) 


where,  as  before: 


f  j(i(M),*---)--s-/T-1-i----r(TrrroxpI"(1/2)£jT(ti,Aj"1<ti)£j<ti)1 


Aj(ti)  -  HjPj(ti-)HjT  ♦  Rj 


— j  ( t  i )  =  ii"  Hjij^i  > 


is  the  dimension  of  z  (number  of  measurements) 


R.  is  the  measurement  noise  strength 


ii 
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II. 3. 2  Sliding  the  Moving  Bank.  When  the  true 
parameter  vector  is  associated  with  a  point  within  the 
moving  bank,  the  moving  bank  operates  esentially  like  a 
smaller  version  of  the  full-bank  estimator.  However,  when 
the  true  parameter  point  lies  outside  the  moving  bank,  this 
condition  must  be  detected  and  some  action  must  be  taken  to 
bring  it  within  the  moving  bank  or  the  estimate  will  be 
erroneous  [28:22,29:10].  Since  the  true  parameter  point  is 
both  unknown  and  uncertain,  some  means  detecting  when  it  is 
not  within  the  moving  bank  must  be  determined.  Addition¬ 
ally,  the  estimator  can  be  expected  to  operate  best  when  the 
moving  bank  is  centered  on  the  true  parameter  point, 
therefore  even  if  the  true  parameter  point  is  within  the 
bank  but  is  close  to  the  perimeter,  the  bank  should  be 
moved.  Maybeck  and  Hentz  [28:22-24,29:10:12]  investigated 
four  means  of  detecting  when  the  true  parameter  point  was 
not  within  the  moving  bank. 

II. 3. 2.1  Residual  MfiSA.t.£>.£jLfiA»  Let  a  likelihood 
quotient  Lj(t^)  be  defined  of  the  quadratic  form 
appearing  in  Equation  (II-5): 

LjUi)  =  rjT(ti)Aj'1(ti)  rjU^  (11-26) 

When  the  true  parameter  point  is  outside  the  moving  bank, 
all  of  the  likelihood  quotients  for  the  elemental  filters 
within  the  moving  bank  can  be  expected  to  exeed  some 
threshold  level  T,  the  numerical  value  of  which  is  set  in  an 


ad  hoc  manner  daring  performance  evaluations.  This 
detection  method  would  indicate  that  the  moving  bank  should 
be  moved  at  time  t^  if: 

min{L1(ti),  LjUi).  ....  LjCt^}  1  T  (11-27) 

The  bank  should  be  moved  in  the  direction  of  the  filter  with 
the  smallest  L-,  as  that  filter  can  be  expected  to  be 
nearest  to  the  true  parameter  point.  This  method  of 
detection  responds  effectively  and  quickly  to  a  real  need  to 
move  the  bankbut  is  also  apt  to  respond  erroneously  to  a 
single  instance  of  large  measnrement  corruption  noise. 

1 1. 3. 2. 2.  Parameter  Position  Est imate  Monitor ina. 

The  current  estimate  of  the  parameter  vector  is  adapted  from 
Equat ion  ( 1 1 —  8 )  : 


J 

A  \ 

•  (tj)  =  L  a  -p  jUj)  (11-28) 

j  =  l  J 

Using  this  logic  would  require  a  move  anytime  the  bank  was 
not  centered  on  the  point  closest  to  the  estimated  parameter 

A 

point.  Since  a.(t^)  depends  on  a  history  of  measurements, 
this  method  of  detection  is  not  as  sensitive  to  a  single 
instance  of  large  measurement  noise  as  is  residual 


moni tor ing . 


11.3.2.3.  Parameter  Position  and  Velocity  Estimate 
Monitor ing.  The  history  of  parameter  position  estimates 
can  be  used  to  generate  an  estimate  of  the  velocity  of  the 
true  parameter  point  in  the  case  where  the  parameters  are 
varying  slowly  and  steadily  in  time.  In  this  case,  the  bank 
is  moved  after  the  current  update  but  before  the  next 
propagation  cycle,  to  center  it  upon  the  point  closest  to 
the  projected  parameter  point.  In  Maybeck  and  Hentz's 
investigations,  this  method  of  detection  was  found  to  perform 
worse  than  did  parameter  position  estimate  monitoring  or 
probability  monitoring  described  below  [28:85,29:23]. 

II. 3. 2. 4  Prob  ability  Monitor ing.  The  conditional 

hypothesis  probabilities  Pj(tj)  computed  via  Equation 
(11-25)  are  monitored,  and  if  the  largest  PjCtj)  is  larger 
than  a  chosen  threshold,  the  bank  is  centered  on  that  filter. 
Maybeck  and  Hentz  [28:85,29:23]  found  that  this  method  of 
detection,  when  used  by  itself,  provided  performance  as  good 
as  parameter  position  monitoring  and  required  less 
comput  a  t ion  . 


1 1. 3. 3.  Changing  .the.  Di^c  r  e.ilia  t„ion  The  filters  in 
the  moving  bank  need  not  be  those  associated  with  adjacent 
points  in  the  parameter  space.  As  seen  in  Figure  II-2,  it 
may  be  more  appropriate  to  space  the  filters  implemented  in 
the  moving  bank  widely  over  the  parameter  space.  This  can 
be  expected  to  decrease  the  accuracy  of  the  estimate,  but 
the  probability  that  the  true  parameter  point  will  lie 


- .  Y.  .-; 
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the  parameter  estimate: 


E{[a  -  £( t i ) ] [a  -  i(ti)]T|z(ti)  =  Z £]  = 

J 

^  [a  -  i(t.)][4j-  a(ti)]Tpj(ti)  (11-29) 
j=l 

drops  below  some  selected  threshold.  The  contraction  from 
the  coarsest  discretization  to  the  finest  discretization  can 
be  done  either  in  one  step  or  several. 

The  bank  can  also  be  expanded  during  normal  operation. 
This  would  be  appropriate  if  the  true  parameter  point 
underwent  a  large  jump  change  in  position  (due  perhaps  to 
some  large  and  sudden  physical  change  in  the  system).  Such 
a  condition  can  be  detected  by  the  residual  monitoring 
discussed  earlier  [28:84,29:20];  in  this  case  all  of  the 
Lj's  would  be  expected  to  be  large  and  close  to  each  other 
in  value.  Probability  monitoring  may  also  give  some 
indication  of  this  condition  by  the  pj's  all  becoming  close 
to  each  other  in  value.  After  the  bank  is  expanded,  it  is 
then  allowed  to  contract  subsequently  about  the  new 
parameter  point,  in  the  manner  discussed  earlier. 

1 1. 3. 4  Ini.t.ia  1. iiAiion.  When  the  moving  bank  is 

moved,  expanded  or  contracted,  it  is  necessary  to 
initialize  any  filters  that  were  not  implemented  in  the 
moving  bank  before  the  action  took  place.  Each  new  filter 
must  be  assigned  values  for  ♦ j ,  B^j,  Kj#  H ^ .  Aj(ti). 
and  Pjft^  [28:28,29:13  ].  All  of  these  except  £j  and  pj 
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bank,  it  would  either  be  three  or  five  filters.  Good 


performance  may  be  achieved  by  dividing  the  total 
probability  weight  of  one  minus  the  sum  of  the  unreset  Pj's 
among  the  new  filters.  This  can  either  be  divided  equally, 
or  it  canbe  apportioned  in  a  manner  indicating  the  new 
filters'  expected  correctness  as  described  in  Equation 
(11-30)  [28:29.29:13]: 


Jc  h 


<»!>  = 


fj(i(ti))(l  -  ^  Pi-Uj)) 
_ finch _ 


(11-30) 


where  ch  =  changed,  unch  »  unchanged,  and: 

fj  (  z(  t  i)  )  = _ 1 _ exp[-(l/2)£jT(ti)Aj-1(ti)ij(ti)] 

(  2  IT )  m/ 2  |  A  j  ( t  i )  I  (1/2) 

Aj(tj),  m,  and  R  are  as  defined  before,  and: 

£j(ti)  =  z.- 


This  may  require  enough  additional  computation  time  that 
there  is  no  net  performance  improvement  over  dividing  the 
net  probability  weight  evenly  among  the  changed  filters. 

For  an  expansion  or  contraction,  it  is  likely  that  all  of  the 
filters  are  changed  or  that  the  old  pj(t{)'s  are  no  longer 
valid,  in  which  case  setting  all  the  Pjftjl's  to  1 / J  is 
appropriate. 


1 1. 3. 5.  Mov  ina-Bank  Adapt ive  Control.  The  controller 

for  nse  with  the  moving-bank  Multiple  model  adaptive 
estimator  can  be  designed  following  any  of  the  three  methods 
discussed  for  the  full-bank  estimator  in  Section  II. 2. 3. 

The  only  difference  is  that  for  the  moving  bank  case,  the 
multiple  model  adaptive  controller  only  uses  the  controllers 
associated  with  the  filters  currently  implemented  in  the 
moving  bank.  Naybeck  and  Bentz  received  good  and 
essentially  identical  results  with  a  single  change ab 1 e-g a  in 
controller  and  a  moving-bank  multiple  model  adaptive 
controller.  The  single  fixed-gain  controller  performed 
poorly  when  the  true  parameter  vector  differed  significantly 
from  the  nominal  parameter  vector  for  which  the  controller 
was  designed.  In  addition,  performance  was  sensitive  not 
only  to  the  magnitude  of  the  error  made  in  estimating  the 
parameter  vector,  but  also  was  more  sensitive  to  errors  in 
one  parameter  than  the  other,  and  whether  the  parameter  was 
overestimated  or  underestimated  [28:104-105,  29:25-27]. 

Naybeck  and  Hentz  also  found  it  necessary  to  shut  off 
control  for  the  initial  period  during  which  the  moving  bank 
was  acquiring  the  parameter  estimate.  If  control  was 
applied  before  the  parameters  were  identified,  the  wrong 
control  was  often  applied,  driving  the  system  unstable 
[28:104,  29:25].  Turning  off  the  control  was  also  necessary 

for  the  same  reason  after  a  jump  change  in  the  true 
parameter  point.  This  led  to  the  recommendation  that  an 
appropriate  method  of  deciding  when  to  enable  control  be 


investigated,  based  on  a  determination  that  parameter 
acquisition  had  taken  place.  The  possibility  of  using  the 
fixed  gain  controller  based  on  a  nominal  parameter  vector 
during  the  acquisition  phase,  rather  than  simply  disabling 
control  entirely  was  also  suggested  [28:106]. 

1 1 . 4  Amfeiiuit.y  Functions  Analysis 

A  tool  that  can  be  used  to  predict  the  ability  of  the 
moving-bank  multiple  model  adaptive  estimator  to  center 
itself  correctly  on  the  filter  in  the  bank  which  is  closest 
to  the  true  parameter  point,  is  ambiguity  function  analysis 
[5:97].  The  generalized  ambiguity  function  is  qiven  by: 


Ai(a.At) 


where  a  is  the  parameter  vector,  a.t  is  the  true  parameter 
vector,  and  L[&,Z^]  is  a  likelihood  function  [5:97-99].  For 
a  qiven  value  of  a.(>  this  function  of  a  yields  information 
about  the  expected  ability  of  the  adaptive  filter  to 
estimate  parameters.  When  plotted  on  a  three  dimensional 
surface,  over  the  plane  of  the  two-dimensional  parameter 
space,  the  curvature  of  the  ambiguity  function  at  At 
predicts  the  precision  with  which  the  adaptive  estimator  can 
estimate  the  parameters.  The  moving  bank  will  center  itself 
on  the  highest  local  peak  within  its  area;  if  there  is  no 
local  peak,  the  moving  bank  will  move  uphill  until  it 


encounters  either  a  local  peak  or  the  edge  of  the 


space.  If  there  are  multiple  peaks  in  the  parameter  space, 
the  moving  bank  may  converge  to  incorrect  parameter  valnes, 
depending  on  the  bank  starting  position.  In  addition,  the 
greater  the  curvature  at  the  peak,  the  greater  is  the 
precision  [28:332].  Examination  of  the  ambiguity  functions 
may  be  useful  in  determining  starting  points  for  the  moving 
bank,  appropriate  discretization  levels,  and  contraction  and 
expansion  strategies.  It  is  important  to  note  that  the 
ambiguity  function  varies  with  £t;  thus  a  different  plot  is 
required  for  each  true  parameter  point  of  interest. 

The  curvature  of  the  ambiguity  function  is  inversely 

✓ 

related  to  the  Cramer  Rao  lover  bound  on  the  estimate  error 
covariance  matrix  by  [5:97]: 

E[[a  -  a.tl  [a  -  AtlT]  l  [-<  dV  3  )  A^,^)  |  ^  (11-32) 

The  ambiguity  function  (Equation  (II-3D)  can  be 
calculated  through  the  evaluation  of  covariance  analysis 
results  in  which  the  true  system  is  based  on  at  and  the 
estimator  is  a  single  Kalman  filter  of  the  same  structure  as 
the  true  system,  but  based  on  a.  instead  of  a.t»  The 
ambiguity  function  can  thus  be  written  [5:98,  28:333]: 

A i ( a , a t )  =  m/2  ln(2n)  -  1/2  In [ I  A ( t £ ; a ) I ] 

-  1/2  tr{A_1(ti;i) [e(ti)Pe(ti‘;at.A)eT(ti)+R(ti)]} 

-  n/2  1  n  (  2  jt  )  -  1/2  In  [  I  P ( t  £  + ; &)  I  1 

-  1/2  tr{P‘1(ti+;i>Pe(ti+;at, a)}  (11-33) 
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where 

ACt.jj.)  =  [H(ti)P(ti~;A)HT(ti)  +  Ktj)]"1 

for  the  Kalman  filter  based  on  & 

+ 

Pe(ti~;it*A^  is  the  covariance  matrix  of  the 

error  between  the  state  estimates  of 
the  Kalman  filter  based  on  &  and  the 
states  of  the  true  system  based  on 
a^,  where  or  '  +  '  denotes  before 

of  after  incorporation  of  the  itb 
measurement. 

'm'  is  the  number  of  measurements, 
and  *n'  is  the  number  of  states. 

The  first  three  terms  are  in  actuality  summed  over  the  last 
N  sample  times;  here  N  is  set  equal  to  one,  which  reduces 
the  size  of  the  fluctuations  in  the  value  of  A^^.a^) 
(flattens  the  plotted  surface).  This  does  not  aid  the 
analysis  of  the  plots,  but  does  make  it  computationally 
s imp  1 e  r . 

I I . 5  Summary 

The  algorithms  for  full-scale  MMAE  and  the  moving-bank 
multiple  model  adaptive  estimator  have  been  developed  in 
this  chapter.  Both  estimators  are  expected  to  give  accurate 
adaptive  estimates  of  true  system  states  in  many  appli¬ 
cations;  however,  the  moving-bank  estimator  is  expected  to 
be  a  more  practical  estimator  for  implementation  because  of 
the  reduced  computational  loading. 


Several  areas  specific  to  the  moving-bank  estimator 
were  explored.  Four  means  of  detecting  when  movement  of  the 
moving  bank  is  required  were  discussed: 

a.  residual  monitoring 

b.  parameter  position  estimate  monitoring 

c.  parameter  position  and  velocity  estimate 
mon i to  r ing 

d.  probability  monitoring 

Changing  the  discretization  of  the  moving  bank  was 
discussed  both  as  an  approach  to  initial  acquisition  of  the 
unknown  parameter  vector  and  as  means  for  r e ac qu i s i t ion 
after  a  jump  change  in  the  unknown  parameter  vector. 
Parameter  estimate  covariance  monitoring  and  probability 
monitoring  were  discussed  as  means  of  detecting  when  the 
bank  should  be  contracted.  Residual  monitoring  and 
probability  monitoring  were  discussed  as  means  of  detecting 
when  the  bank  should  be  expanded. 

When  moving,  expanding,  or  contracting  the  moving  bank, 
it  is  necessary  to  reset  the  state  estimates  and  probability 
weightings  of  the  new  filters  in  the  bank.  The  current 
overall  state  estimate  is  appropriate  for  resetting  the 
individual  state  estimates.  After  the  bank  has  been 
expanded  or  contracted,  it  is  appropriate  to  reset  the 
probability  weightings  of  all  of  the  filters  to  1/J.  After 
the  bank  has  been  moved,  only  the  probability  weightings  of 
the  new  filters  should  be  reset  by  dividing  among  them  the 
total  probability  weight  left  after  subtracting  the 


probability  weightings  of  the  unchanged  filters  from  the 
total  weight  of  one.  The  remaining  weight  can  be  divided 
either  evenly  or  in  a  manner  reflecting  the  estimated 
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correctness  of  the  new  filter. 

This  research  is  intended  to  explore  and  expand  these 
movement,  expansion,  and  contraction  decision  making  and 
reset  methods,  as  applied  to  the  control  of  the  bending 
modes  in  a  realistic  model  of  a  large  flexible  space 
structure.  The  model  of  the  flexible  space  structure  is 
developed  in  the  following  chapter. 


III. 


SPACE  STRUCTURE  MODEL 


1 1 1 . 1 .  Introduc  t ipn 

This  chapter  describes  the  system  equations  for  the 
'Draper  Labratory  /  Rocket  Propulsion  Labratory  (RPL) 
Configuration'  model  of  a  large  flexible  space  structure, 
illustrated  in  Figure  III— 1.  The  structure  consists  of  a 
rigid  central  hub  with  four  structurally  identical  flexible 
appendages  cantilevered  radially  from  the  hub.  Experimental 
work  concerning  the  planar  rotational/  vibrational  dynamics 
of  this  U.S.A.F.  Rocket  Propulsion  Laboratory  demonstration 
model  is  being  conducted  at  the  Charles  Stark  Draper 
Laboratory,  where  a  physical  model  of  the  structure  is 
supported  on  an  air  bearing  table  [31], 

The  differential  equations  describing  the  model  are 
developed  from  the  equations  of  motion  for  the  unforced 
system  by  Muckenthaler  [32];  the  expressions  for  the  kinetic 
and  potential  energy  are  developed  and  discretized,  using 
the  assumed  modes  method  with  terms  of  higher  order  than  two 
ignored.  This  allows  the  mass  and  stiffness  matrices  to  be 
identified.  The  integral  equations  for  the  mass  and 
stiffness  matrices  are  solved  symbolically  and  a  closed  form 
expression  for  them  is  achieved.  The  second  order 
differential  equations  are  then  placed  in  the  standard  state 
space  form  for  a  stochastic  system  model.  This  allows  the 
construction  of  the  elemental  Kalman  filters  used  in  this 
thesis  for  the  investigation  of  moving-bank  MMAE. 
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Figure  III-l.  Draper/RPL  Configuration  Model 
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1 1 1 . 2  .  Second  Or  4q  £  Form 

As  described  by  Hackenthaler  [32 ]  the  motion  of  the 
uncontrolled  Draper/RPL  configuration  model  can  be  described 
by  the  second  order  vector  differential  equation  [32:25] 

Mq+Cq  +  Kq=0  (1II-1) 


where  : 

H:  is  the  N  x  N  symmetric  and  positive  definite  mass 

matr  ix 

C:  is  the  N  x  N  structural  damping  matrix 

K:  is  the  N  x  N  symmetric  positive  semi-definite 

stiffness  matrix 

q:  is  the  N  dimension  vector  of  generalized  coordinates 
q  =  [  e .  oa .  .  .  Op  ,  vx .  .  .  Vp  ]  T 

where  6  is  the  rotational  displacement  of  the  hub. 
and  Uj  and  Vj  are  the  generalized  displacement 
coordinates  of  the  arms,  and  N  =  2p  +  1,  where  p  is 
the  number  of  modes  of  interest. 

1 1 1. 2.1.  Aqqqmpt iflSA.  The  equations  developed  to 

describe  the  Draper/RPL  configuration  model  are  based  on  the 
following  assumptions  [32:11]: 

a)  The  longitudinal  and  out-of-plane  vibrations  of  the 
arms  are  of  much  higher  frequency  than  the  transverse 
vibrations  and  are  negligible. 

b)  An t i- s y m m e t r  i  c  deformations  are  such  that  the 
deflections  of  the  first  arm  are  equal  in  magnitude  but 


opposite  in  direction  to  the  second  arm,  and  the  deflections 


Figure  III-2.  Anti-Symmetric  Deflections  [32:11]. 

of  the  third  arm  are  equal  and  in  magnitude  but  opposite  in 
direction  to  the  fourth  arm;  see  Figure  IXI-2. 

c)  Since  the  hub  of  the  physical  model  is  supported  by 
an  air  bearing  table,  the  force  of  gravity  can  be  neglected 
in  the  derivation  of  potential  energy. 

d)  The  structural  damping  coefficients  are  very  small; 
therefore  the  structural  damping  matrix  C  is  approximated  as 
the  zero  matrix. 

e)  All  four  arms  are  structurally  identical  with: 

Lj  —  L  2  =  L  3  =  L 

*1  =  l2  =  fc3  =  *4  =  1 

h  ^  =  1*2  =  h  ^  =  ^4  22  ^ 

and  the  four  tip  masses  are  identical  and  equal  to  m2* 


III. 2. 2. 


M and  S  t.  Eanat  ions.  The  mass  and 
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stiffness  matrices  can  found  by  solving  for  the  kinetic  and 
potential  energy  of  the  free  vibation  system.  Neglecting 
the  terms  of  order  three  or  higher,  the  kinetic  and  potential 
energies  can  be  written  in  the  form  [32:23]: 

T  =  1/2  flT  M  i 
V  =  1/2  aT  K  a 

where  the  elements  of  the  mass  matrix  are  given  for 
i  =  l,...,p  ;  and  j  =  l,...,p  ;  by: 

M(l.l)  =  2  (Ic  +  I2<J) 

R+L 

M ( 1 , i  +  1 )  =  2 p th  I  x0£  dx  +  2m2(R+L)0i( z=x-R=L) 

R 

+  2 I2c0 ’ j ( z=x-R=L) 

R+L 

M(I,i+p+l)  =  2pth  J  y0 j  dy  +  2m2 ( R+L) 0 ^ ( z  =  y-R=L) 

R 

+  2I2(j0  '  i(z  =  y-R=L) 

M (  i  +  1,1)  =  M(l,  i  +  1) 

M (  i  +  p  +  1 , 1 )  =  M ( 1 ,  i  +  p  +  1 ) 

R+L 

M ( i  +  1, j+1)  =  2  p  t  h  j  0i0j  dx 

R 

+  2mo0 ; ( z=x-R=L)0 ; ( z=x-R=L) 

+  1  J 


+  2I,„  0'1(z  =  x-R=L)0*  .(z  =  x-R=L) 


M ( i  +  p  +  1 , j +p  +  l )  =  2  p  t  h  J 


0i0 j  dy 


R 

+  2a»20  A  (  z  =  y~R=L)  0  j  (  z  =  y-R=L) 

+  2I2c0'i(z  =  y-R=L)0' j  (  z  =  y-R=L) 

M (  i+1 , j  +p  +  l )  =  0 

M ( i+l+p,  i  +  1)  =  0 

where  0^  is  a  function  describing  the  bending  modes  of  the 
arms,  defined  as: 

0^2)  =  1  -  cos  (  inz/L)  +  (1/2)  (-1)  i  +  1(  inz/L)2 
z=x-R(ory-R) 

and  01 j  is  the  first  derivative  of  0j  with  respect  to  z. 

The  pysical  dimensions  of  the  model  t,  h.  m2,  I2c*  1 0 . 
R.  and  L,  are  as  shown  in  Fignre  III-l.  In  addition  p  is 
the  mass  density  of  the  arms,  and: 

Ic  =  ( 1 / 2 ) I Q  +  2I1  +  2m2 ( R+L) 2 

R+L 

Ij  =  pth  J  x2  dx 
R 

The  elements  of  the  stiffness  matrix  are  given  for 
i  =  l,...,p  ;  and  j  =  l,...,p  i  by: 

1(1.1)  =  0 
K(  i  +  1 , 1 )  =  0 
K (  i  +  l+p , 1 )  =  0 


K  ( 1 ,  j +p+l )  =  0 


L 

K(  i  +  1 , j  +  1)  =  2EI  J  0i"(z)0j"(z)  dz 

0 

K (  i+l+p , j +1 +p )  =  K(i+l,j+l) 

K(i+l+p,j+l)  =  0 

K (  i  +  1 ,  j  +l+p )  =  0 

where  E  is  the  M  o  d  u  1  u  s  of  Elasticity  of  the  arms,  and  I  is 
the  area  moment  of  inertia  based  only  on  t  and  h  as: 

I  =  (l/12)ht3 

and  0"jL  tlie  second  derivative  of  0^^  with  respect  to  z. 

The  equations  for  the  elements  of  the  N  and  K  matrices 
are  simplified  and  evaluated  in  Appendix  A. 

1 1 1. 2. 3.  Eigenvalues  of  the  free  system.  The 

eigenvalues  of  the  free  system  can  be  found  by  solving  the 
generalized  eigenvalue  equation: 

Kx  =  IMx  ( I I 1-2 ) 

The  square  roots  of  the  eigenvalues  represent  the 
vibrational  frequencies  of  the  modes  depicting  the  relative 
motion  of  the  free  arms.  The  modes  appear  in  pairs,  very 
close  in  frequency,  depicting  what  Junkins  [31]  calls 
'unison'  and  'opposition*  modes.  The  opposition  modes  are 
simple  cantilever  beam  modes  characterized  by  the  adjacent 
beams  moving  in  opposition.  He  states  the  unison  modes  are 
perturbed  cantilever  modes,  with  just  slightly  higher 
frequencies,  with  all  four  beams  moving  in  unison  and  the 
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Opposition  Mode 


Uni  son  Mode 


Figure  III-3.  Opposition  and  Dnison  Modes  [32:26]. 

hub  having  non-zero  rotation.  The  hub  rotation  is  necessary 
to  counter  the  unison  movement  of  the  arms  and  tip  masses 
to  conserve  angular  momentum  in  the  system  [32:26];  see 
Figure  II 1-3. 

1 1 1 . 3  Sjt^t  e  Sj>  a£ e.  Farm 

The  second  order  differential  equation  describing  the 
free  system.  Equation  (III-l),  can  be  modified  to  describe 
the  controlled  system  as: 

Ma+Ka=Bu  ( III-3 ) 

where  B  is  anN  x  M  control  matrix  and  a  is  an  M-dimensional 
control  vector  added  to  the  original  Equation  (III-l). 

Recall  that  structural  damping  was  assumed  to  be  very  small 
so  the  damping  matrix  C  has  been  neglected. 
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then  Equation  ( 1 1 1  —  3 )  can  be  rewritten  as: 


-  M~ 1 K  0 


in  first  order  fora, 


-M-1B 


( I I 1-4 ) 


On  the  Draper/RPL  configuration  model,  actuators  are 
located  one  on  the  hub  and  one  halfway  along  each  arm. 

These  apply  a  torque  Uj  to  the  hub,  a  torque  U2  at  z  =  L/2 
to  arms  1  and  2,  and  a  torque  u^  at  z=L/2  to  arms  3  and  4. 
Recall  Figure  III-2,  in  which  the  arms  are  numbered.  Matrix 
B  then  becomes  [32:30]: 


B  =  0  2ft' A( z=L/2)  0 

.  0  0  2ft'  i(z  =  L/2)_ 

where  ft^(z)  is  of  dimension  p  (the  number  of  modes  of 
interest).  B  is  therefore  an  N  x  3  matrix. 

Position  and  velocity  measurements  of  the  model  are 
available  from  colocated  position  and  velocity  sensors  on 
the  hub  and  at  positions  along  each  arm.  For  the  state 
space  model,  the  measurement  equation  (assuming  for  the 
moment  noise-free  measurements)  is  thus: 
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£  =  ft 

0  B 
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where  the  measurement  matrix  fl  is  [32:30] 
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1  column  p  columns  p  columns 
H  is  thus  a  5  x  N  matrix. 

1 1 1 . 4  .  Ffirm 

The  mathematical  model  can  be  placed  in  stochastic  form 
by  adding  a  noise  matrix  Gv>  multiplying  the  control  input 
noise  w,  to  Equation  (III-4),  and  a  measurement  noise  vector 
to  Equation  (III-5).  If  it  is  assnmmed  that  input  noise 
will  enter  through  the  actuators  and  measurement  noise 
through  the  sensors,  then  Gw  becomes  identical  to  the 
augmented  B  matrix.  The  system  is  thus  described  by: 
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I I I . 5  Physical  Syilflm  CflflAia fli* 

The  nominal  constants  describing  the  physical  system 
are  listed  in  Table  III-l.  For  the  purpose  of  setting  up  a 


m 


parameter  space  in  which  to  test  moving-bank  MMAE 
algorithms,  p,  the  mass  density  of  the  arms,  and  E,  the 
modulus  of  elasticity  for  the  arms,  were  allowed  to  vary  - 
20%  to  +16%  from  their  nominal  values  in  discrete  steps  of 
4%,  yeilding  a  10  by  10  point  parameter  space.  This  is  an 
unre a  1 i s t ic ly  large  amount  of  variation  in  these  parameters, 
but  it  proved  to  be  necessary  for  the  purpose  of 
investigating  moving-bank  MMAE  algorithms  as  applied  to  this 
system. 


Table  III-l.  Configuration  Constants  for  Draper/RPL  Model 


R: 

hub  radius 

* 

1  ft 

hub  rotary  inertia 

8  slug-ft2 

p: 

mass  density  of  the  arms 

= 

5.22  slug-ft2 

E: 

arm  modulus  of  elasticity 

= 

1.584E09  lb/ft2 

i 

t : 

arm  thickness 

= 

0.0104166  ft 

h: 

arm  height 

= 

0.5  ft 

i 

L: 

arm  length 

= 

4.0  ft 

i 

m2 

:  tip  mass 

= 

0.156941  slug 

1 

L _ 

^2c"  tip  mass  rotary  inetria 

0.0018  slug-ft2 

III. 

6. 

State  Reduc  t ion 

I  f 

p  is  taken  to  be  two  (two 

oppos it  ion 

and  two  un  i  s o n 

modes) 

then  N  is  five,  and  this  results 

in  a 

ten- state  model 

for 

the 

augmented  system.  A  ten- 

state 

model 

is  needlessly 

more 

complex  than  is  required  for 

this 

the  s  i 

s.  If  the 

hypothetical  purpose  of  this  e s t i ma to r / c on t r o  1 1 e r  is  to 
eliminate  uncommanded  torque  to  the  hub  of  the  spacecraft 
(as  might  be  required  if  it  is  desired  to  point  an  apparatus 
on  the  hub  in  a  specified  direction),  then  it  is  only 


necessary  to  control  those  bending  nodes  which  produce  a 
torque  on  the  hub. 

As  discussed  earlier,  only  the  nniion  nodes  produce  a 
torque  on  the  hub  [32:26].  Junkins  [31]  deternined  that  the 
unison  modes  were  the  even  nunbered  nodes  (i.e.  2,4,6,...). 
If  the  opposition  nodes  are  assunned  to  be  stable  the 
natrices  describing  the  systen  can  be  sinplified  by 
eliminating  those  state  variables  depicting  the  positions 
and  velocities  of  the  opposition  nodes.  The  state  vector  is 
then  reduced  to  six  states  and  the  natrices  describing  the 
system  need  to  be  simplified  by  eliminating  the  rows  and 
columns  pertaining  to  the  opposition  modes.  Thus  for  the 
F  matrix,  every  other  row  and  column  is  eliminated,  and 
the  result  is  a  six  by  six  matrix.  For  the  augmented  B  and 
Gw  matrices,  every  other  row  is  eliminated  making  them  six 
by  three  matrices.  For  the  augmented  H  matrix,  every  other 
column  is  eliminated,  leaving  a  ten  by  six  matrix.  Sample 
matrices  for  a  nominal  case  are  shown  in  Appendix  A. 

I I I . 7 .  Summary 

This  chapter  developed  equations  describing  a  realistic 
model  of  a  large  flexible  space  structure.  The  stochastic 
state  space  mathematical  model  derived  is  dependent  on  the 
physical  parameters  describing  the  real-world  system. 
Variation  in  the  physical  parameters  will  create  variation 
in  the  stochastic  model,  and  this  will  allow  investigation 
into  the  use  of  moving-bank  NMAE  to  estimate  both  the  system 
states  and  the  varying  parameters  of  the  physical  system. 
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IV.  SiBBllliQfl 


la  order  to  tTilsite  the  performance  of  the  moving-bank 
multiple  model  adaptive  e s t i ma to r / c on t ro 1 1 e r ,  it  is 
necessary  to  simulate  the  space  structure  and  the 
e s t ima t o r / cont r o 1 1 e r  operation.  The  computer  simulation 
used  provides  a  Monte  Carlo  analysis  of  both  the  space 
structure  model  and  the  e s t i ma to r / c on t r o 1 1 e r.  This  chapter 
contains  a  brief  discussion  of  Monte  Carlo  analysis,  an 
outline  of  the  computer  software  used  to  accomplish  the 
analysis,  and  the  plan  for  analyzing  the  performance  of  the 
moving-bank  algorithm  and  the  specific  logics  used  for  the 
move,  contract,  and  expand  decisions. 

IV. 2.  Mon!*  Cgilfi  Aafilxfiii 

The  random  nature  of  the  input  and  measurement  noise 
processes  makes  it  impossible  to  select  a  single  typical 
example  of  them.  Thus  in  order  to  characterize  the 
performance  of  the  moving-bank  multiple  model  adaptive 
e s t i ma to r / c on t r o  1 1 e r  statistically,  it  is  necessary  to 
examine  the  ensemble  average  of  the  e s t i m a t o r / c on t r o 1 1 e r 
performance  using  many  samples  of  the  error  process.  Monte 
Carlo  analysis  does  precisely  this:  a  number  of  individual 
time  simulations  of  the  e s t i ma t o r / con t r o  1 1 e r  are  made  and 
sample  statistics  (means  and  covariances)  are  computed 
directly  for  each  sample  time  [33:29]. 
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For  the  analysis  performed  in  this  thesis,  the  true 


system  model  is  described  by  the  linear  t ime- inv ar ian t 
stochastic  difference  equation: 


iUj)  =  +  Bda(  t  +  G^U^)  (IV-1) 


where 

b  is  the  state  transition  matrix  from  t^_^  to  t^ 

Bd  is  the  control  input  matrix 
Gd  is  the  noise  input  matrix 

Recall  that,  for  the  Draper/RPL  model  as  shown  in  Equation 
(II1-6),  the  noise  input  matrix  is  identical  to  the  control 
input  matrix,  therefore: 


Bd  =  Gd 

where  Bd  is  the  discrete-t ime  equivalent  of  the  augmented  B 
matrix  of  Equation  (III-6),  given  by  [33:171]: 


Bd  =  j  #(t£,T)B  dx 

t  i_j 

and  £(t^)  is  the  discrete  time  equivalent  of  the  augmented 
state  vector  Q  of  Equation  (III-6).  Also  note  that  (>, 
and  Bd  (and  therefore  Gd)  are  all  functions  of  the  true 
parameter  vector  a.t,  where 


^t  = 


P 

E 


t 

t  - 


The  model  is  driven  by  both  the  known  control  q(t^_^)  and  a 
discrete  zero-mean  white  Gaussian  noise  w(t^_j)  of 
covariance  Qd.  Noise-corrupted  measurements  of  the  system 


states  are  provided  to  tlie  estimator  in  the  form  of: 


ift^  =  +  v(t£> 


( I V-2 ) 


where  H  is  now  the  augmented  version  of  the  measurement 
matrix  shown  in  Equation  (III-7)  and  the  measurements  are 
corrupted  by  a  discrete-time  zero-mean  white  Gaussian 
measurement  noise  y(t^)  of  covariance  R. 

The  true  system  and  the  e s t i m a t or / con t r o 1 1 er  are 
operated  from  time  t0  to  time  tf  for  a  sufficient  number  of 
runs  that  the  computed  sample  means  and  covariances  of  the 
random  variables  of  interest  are  good  approximations  to 
ensemble  averages  (expectations).  It  is  possible  to 
determine  the  number  of  runs  that  is  sufficient  by  observing 
how  the  computed  sample  statistics  change  as  the  number  of 
runs  is  increased;  after  a  sufficient  number  of  runs,  the 
sample  statistics  will  converge  to  a  constant  [33:29].  For 
this  problem  that  entails  on  the  order  of  10  or  more  runs. 

Figure  IV-1  depicts  the  simulation  of  the  true  system, 
the  estimator,  and  the  controller.  The  variables  of 
interest  are: 


the  system  ('truth  model')  states  -  x.t(ti) 


the  error  in  the  estimate  of  the  system  states  - 


S.x<  t  A)  =  it  ( 1 1)  -  x  ( t  £) 


msm 


truth 

system 


i<ti) 

*  i  M  esti-  — 


Figure  IV-1.  System,  Estimator,  and  Controller  Simulation. 


the  error  in  the  parameter  estimate 
£a(ti>  *  *t<*i>  ’ 


and  the  control  input  -  u(t^) 


The  means  and  covariances  of  the  variables  of  interest  are 
all  computed  similarly;  for  example,  the  mean  of  the  error 
in  the  state  estimate  is  computed  as: 


^ t ( t i ) )  r  Mex(ti>  =  (1'N>  ^  &xk(ti) 

k  =  l 


(IV-1) 


where  N  is  the  total  number  of  runs  performed  in  the  Monte 
Carlo  analysis,  and  £.xk^*i^  is  the  value  of  £x(t£)  during  the 


run. 


The  covariance  of  4.x(t£>  is  computed  as: 


Ede^tj)  -  E{£x(ti)}  Hfix(ti)  -  E{ex(ti)}]X)  r  pex(ti)  = 


\  T 

( 1  /  ( 1~N)  i  exk<ti)exk(ti) 


(1/1-N))  Mex<ti)Mex(ti)  ( I V-4  ) 


When  evaluating  the  estimator  alone  the  controller 
block  in  Figure  IV-1  is  replaced  by  a  dither  signal  which 
excites  the  system  states,  bat  is  independent  of  the 
performance  of  the  estimator.  The  dither  signal  is  used 
based  on  the  experiance  of  Hentz  [28:58],  who  found  that 
system  iden t i f iab i 1 i ty  was  enhanced  by  sufficiently  and 
persistently  exciting  the  true  system  modes  with  a  known 
periodic  input.  Without  the  dither,  the  estimator  had 
considerable  difficulty  in  identifying  the  uncertain 
parameters,  and  little  basis  for  evaluating  various  decision 
algorithms  was  found  because  the  system  would  essentially 
reach  a  quiescent  state  reqardless  of  the  algorithm  used 
[28:57-58].  The  specific  dither  signal  to  be  used 
(magnitude,  frequency,  etc.)  is  determined  experimentally  as 
a  rigorous  determination  of  the  optimal  input  to  enhance 
system  id e n t  i  f  i  ab  i  1  i  t y  is  beyond  the  scope  of  this  effort. 

It  is  most  appropriate  to  look  at  the  statistics  of  the 
error  in  the  state  estimate  and  the  error  in  the  parameter 
estimate  in  evaluating  the  performace  of  the  estimator.  The 
error  in  the  state  estimate  gives  the  best  means  of 


comparing  the  estimator  to  other  types  of  estimators,  and 
the  primary  reason  for  adaptation  is  to  enhance  state 
estimation  precision,  rather  than  to  provide  accnrate 
parameter  estimates  for  their  own  sake.  The  characteristics 
of  the  error  in  the  parameter  estimate  reveals  the  accuracy 
of  the  parameter  estimates  that  may  later  be  fed  to  the 
controller.  The  error  in  the  parameter  estimate  can  also 
give  insight  into  how  good  a  job  of  identifying  the  closest 
elemental  filter  the  moving-bank  is  doing,  thereby  giving  a 
means  of  evaluating  the  various  move,  contract,  and  expand 
algorithms;  statistics  on  the  location  of  the  center  of  the 
bank  are  also  valuable  in  this  evaluation. 

When  evaluating  the  e s t i m a t o r / c on t r o 1 1 e r  combination  it 
is  more  appropriate  to  look  at  the  statistics  of  the  true 
state  values  and  commanded  controls.  In  this  thesis,  the 
object  of  the  control  input  will  be  to  drive  the  states  to 
the  quiescent  state;  therefore,  deviations  from  zero  are 
undesirable  characteristics  to  be  analyzed  in  evaluating 
controller  performance.  It  may  also  be  useful  to  examine 
the  control  input  to  detect  unreasonable  commanded  control 
levels. 

IV.  3.  Software  Dej.crifiiifln 

The  Monte  Carlo  analysis  of  the  moving-bank  estimator/ 
controller  required  the  development  of  three  computer 
programs  (for  a  detailed  discussion,  see  Appendix  C.).  The 
first  program  is  a  preprocessor  which  creates  the  parameter 


space.  The  second  program  is  a  primary  processor  which 
perfoms  the  Monte  Carlo  simulation  runs  and  generates  the 
data  describing  each  rnn.  The  third  program  is  a 
postprocessor  which  compotes  the  means  and  covariances  of 
the  variables  of  interest  and  generates  plots  based  on  their 
values.  In  addition,  a  program  which  computes  the  ambiguity 
functions  and  generates  their  plots  was  developed  as  a 
separate  analytical  tool,*  it  requires  a  variation  of  the 
preprocessor  to  set  up  the  parameter  space  values  it  uses. 

The  preprocessor  sets  up  the  parameter  space.  That  is, 
for  each  realization  of  the  uncertain  parameter  vector,  it 
computes  the  matrices  necessary  to  describe  that  parameter 
point  uniquely  in  the  equations  used  to  describe  the  truth 
model  and  filters.  This  requires  that,  for  the  kth 
parameter  point,  the  $k.  Bdk,  Hk,  Kk»  Pj^tj-),  Ak_1,  and 
the  Gk*  matrices  and  the  determinant  of  Ak  be  computed.  The 
preprocessor  takes  as  input  the  ten  discrete  values  each 
parameter  can  assume,  and  the  weighting  matrices  Wz,  and  WQ 
which  describe  the  quadratic  cost  function  used  in  the 
design  of  the  LQG  controller  gain  (G*)  matrices.  The 
preprocessor  outputs  a  file  containing  all  the  matrices 
calculated,  the  weighting  matrices,  the  strength  of  the 
input  noise,  and  the  covariance  of  the  measurement  noise. 
Unless  the  inputs  change,  the  file  containing  the  output  can 
be  used  for  all  of  the  simulation  runs,  making  it 
unnecessary  to  rerun  the  preprocessor. 


The  primary  processor  performs  the  Monte  Carlo 
simulations.  It  consists  of  an  executive  routine  which 
calls  several  subroutines.  After  input  and  initialization, 
the  executive  routine  consists  of  an  outer  loop  that  counts 
the  desired  number  of  Monte  Carlo  runs,  and  an  inner  loop 
that  performs  the  operations  necessary  for  each  sample 
period  from  the  starting  time  of  the  simulation  (tQ)  to  the 
ending  time  (tf).  For  each  sample  period,  separate  routines 
are  called  to  propagate  the  true  system  from  the  last 
sample  time,  propagate  the  filters  currently  implemented  in 
the  moving  bank  from  the  last  sample  time,  take  a  noise- 
corrupted  measurement  of  the  true  system,  update  the  filters 
in  the  moving  bank,  calculate  the  control  input,  and  finally 
make  the  decision  whether  to  move,  expand,  or  contract  the 
bank.  If  it  is  decided  to  alter  the  bank,  separate  routines 
are  called  to  perform  the  move,  expansion,  or  contraction. 
After  the  computation  for  each  sample  period  is  complete, 
the  variables  of  interest  are  written  to  a  data  file. 
Following  each  run,  the  outer  loop  reinitializes  all 
necessary  variables  before  begining  another  run.  Inputs  to 
the  primary  processor  are  the  file  containing  the 
description  of  the  parameter  space,  true  system  parameters, 
a  mode  input  which  specifies  the  move/ contr ac t/expand 
algorithms  to  be  used  and  associated  thresholds,  initial 
probability  weightings  for  the  filters  in  the  moving  bank, 
and  initial  filter  states.  The  output  of  the  primary 
processor  is  a  data  file  for  each  variable  of  interest 


covering  all  of  the  runs,  and  a  more  detailed  print  file 
covering  jnst  the  first  run  which  includes  information 


describing  the  exact  filters  implemented  in  the  moving  bank 
and  the  variable  affecting  the  decision  algorithms.  The 
print  file  allows  detailed  analysis  of  unexpected  or  unusual 
results,  but  avoids  excessive  output. 

The  post-processor  takes  a  data  file  of  a  variable  of 
interest  and  calculates  the  sample  means  and  covariances  for 
each  sample  time  from  tQ  to  tf.  The  post-processor  then 
generates  plots  of  time  histories  of  the  means  of  each 
variable  of  interest  +  la,  where  a  is  its  standard 
deviation.  Thus  each  data  file  generated  by  the  primary 
processor  requires  a  separate  rnn  of  the  post-processor; 
this  provides  for  simplicity  and  flexibility  in  determining 
which  variables  to  plot. 

The  simulation  requires  that  the  driving  noise  wlt^^) 
and  the  measurement  noise  v(t^)  be  zero-mean  white  Gaussian 
processes.  Fortran  provides  a  random  number  generator  which 
can  be  used  to  approximate  the  required  random  variable’s 
realizations  at  each  time.  If  yi  is  the  random  variable 
available  directly  from  the  random  number  generator, 
uniformly  distributed  between  0  and  1,  then  a  zero-mean 
Gaussian  random  variable  with  a  variance  of  1  can  be 


In  order  to  simulate  a  zero-mean  Gaussian  random  vector  with 


covariance  Qd,  the  following  operation  is  performed: 

w  =  Dr 

where  the  elements  of  r.  are  computed  by  independent  calls  to 
Equation  (IV-1)  for  each  scalar  component,  and  where 

D  =  \/Qj;  i.e.  Qd  =  DDT 

For  the  sake  of  simulation,  the  Cholesky  decomposition  is 
used  to  generate  the  square  root  [33:408]. 

IV. 4.  S imul a t ion  Plan 

A  systematic  approach  will  be  used  to  study  the 
performance  of  the  moving-bank  e s t i m a t or / con t r o 1 1 e r.  The 
performance  analysis  will  be  divided  into  two  main  parts. 
First  the  performance  of  the  estimator  alone,  without 
feedback  control  (i.e.,  the  controller  block  in  Figure  IV-1 
is  replaced  with  a  dither  signal  independent  of  the  state 
estimates),  is  evaluated.  Analysis  in  this  portion  of  the 
study  will  concentrate  on  qualitatively  identifying  the 
'best'  estimator  configuration  which  will  then  be  used  for 
the  analysis  of  the  controllers  to  identify  the  'best' 
adaptive  e s t i m a t o r / c on t r o 1 1 e r  combination. 

The  performance  of  the  estimator  alone  is  accomplished 
by  driving  the  true  system  with  a  zero-mean  white  noise  in 
combination  with  a  dither  signal  for  each  of  the 
move / c on t r ac t / e xp and  decision  algorithms  to  be  tested. 

For  each  decision  algorithm  being  evaluated,  the  true 
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parameter  vectors  will  be  chosen  which  exercise 

that  algorithm.  The  estimator  will  be  evaluated  first  using 
only  movement  at  the  finest  discretization  level,  and  the 
issue  of  initializing  the  probability  weightings  of  the  new 
filters  added  to  the  bank  will  be  explored.  Then  the 
effects  of  contraction  from  a  coarse  discretization  to  a 
fine  discretization  will  be  explored  and  contrasted  to 
straight  movement.  If  contraction  proves  beneficial, 
expansion  to  a  coarse  discretization  following  the  detection 
of  a  jump  change  in  the  trio  parameter  point  will  be 
investigated. 

To  evaluate  the  effect  of  computing  probability 
weightings  for  new  filters  added  to  the  bank  during  a  move 
based  on  their  expected  'correctness'  (e.g.,  based  on  how 

far  the  new  filters'  g.  values  are  from  the  current 
aCt^))  vs.  setting  the  probability  weightings  of  the  new 
filter  to  an  equal  share  of  the  total  probability  of  the 
filters  removed,  the  system  uncertain  parameters  will  be 
constant  over  tQ  to  tf,  and  set  equal  to  one  of  the 
discretized  points  in  the  parameter  space  different  from  the 
initial  bank  center.  The  effect  of  the  probability 
weightings  can  then  be  seen  in  the  speed  with  which  the 
parameter  estimate  converges  to  the  true  system  parameter 
point,  and  in  the  error  in  the  state  estimates.  Movement 
primarily  takes  place  when  the  true  parameter  point  lies 
outside  the  area  encompassed  by  the  moving  bank,;  it  does  not 
matter  whether  the  true  parameter  point  is  exactly  equal  to 


been  a  catastrophic  failure  in  one  of  the  aris  of  the  space 
strnctore  model)  than  would  be  possible  by  allowing  the 
estimator  to  converge  to  the  new  point  by  movement  alone. 
Therefore,  the  results  of  the  simulations  using  expansion 
will  be  compared  to  simulations  using  the  same  jump  changes 
in  parameter  value  where  expansion  is  not  allowed. 

After  the  evaluation  of  the  individual  estimator 
algorithms,  a  composite  'best'  estimator  will  be  determined, 
made  up  of  the  move / con t rac t / expand  algorithms  that 
performed  the  'best'.  The  determination  of  'best'  will  be 
based  upon  a  tradeoff  between  added  computational  loading 
and  faster  state  and  parameter  aquisition  times,  lower  state 
and  parameter  estimate  biases,  and  lower  state  and  parameter 
estimate  error  variances.  The  goal  in  moving  from  full- 
scale  NMAE  to  moving-bank  MNAE  is  to  obtain  an  estimator 
with  similar  performance  but  which  has  enough  of  a  reduction 
in  the  computational  loading  required  to  make  it  practical 
for  more  applications.  Therefore,  any  decision  algorithm 
which  increases  computational  loading  is  a  step  in  the  wrong 
direction  unless  this  additional  loading  is  outweighed  by 
significant  gains  in  performance. 

When  the  best  composite  estimator  has  been  determined, 
that  composite  estimator  will  be  used  to  evaluate  the 
controller  configurations  disscussed  in  Chapter  II: 

a.  single  fixed-gain  controller 

b.  single  changeable-gain  controller 

c.  moving-bank  multiple  model  adaptive  controller 


The  «stin*tor/ controller  combination  will  be  evaluated  using 
both  constant  true  system  uncertain  parameter  values  and 
jump  changes  in  the  true  system  uncertain  parameter  values. 
In  addition,  an  evaluation  will  be  made  of  the  effect  of 
turning  off  control  during  the  parameter  acquisition  phase 
vs.  using  the  the  adaptive  controller  from  time  tQ»  or  a 
nominal  fixed  gain  controller  until  parameter  acquisition 
takes  place  then  transferring  to  the  adaptive  controller. 

As  in  the  estimator-only  case,  any  increase  in  computational 
complexity  must  be  offset  by  significant  gains  in 
performance  (lover  state  biases  and  variances  from  the 
quiescent  state)  in  order  to  justify  the  decrease  in 
practicality  associated  with  the  increase  in  computational 
loading . 

IV . 5  Summary 

This  chapter  has  discussed  the  overall  method  of  Monte 
Carlo  analysis,  the  specific  organization  of  the  software 
used  for  the  simulation,  the  plan  and  criteria  for 
evaluation  of  the  various  estimation  and  control  algorithms, 
and  associated  move,  contract,  and  expand  logics.  It  was 
noted  that  increases  in  the  e s t i m a t o r / c on t r o 1 1 e r  performance 
at  the  expense  of  increased  computational  loading  must  be 
significant  to  justify  the  decrease  in  practicality.  The 
results  of  the  simulations  are  presented  in  the  next 
chapter. 
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RESULTS 


V.  I  I,at£fldBgt.lfiR 

This  chapter  presents  the  results  of  the  Monte  Carlo 
simulations.  The  goals  of  the  simulations  were  to  evaluate 
the  effectiveness  of  moving-bank  multiple  model  adaptive 
e s t i ma t ion/ con t ro 1  when  applied  to  the  realistic  situation 
described  in  Chapter  III,  and  specifically  to  evalulate  the 
various  move,  contract,  expand,  and  algorithms  described  in 
Chapter  II.  These  goals  vere  only  partially  met  due  to 
numerical  difficulties  which  forced  the  use  of  approx¬ 
imations  in  computing  probability  weightings,  and  the 
ambiguity  function.  In  addition  the  system  being  controlled 
proved  not  to  require  adaptive  control,  despite  the 
uncertainties  in  the  parameters  describing  it.  This  made 
the  evaluation  of  the  control  algorithms  impossible. 

V . 2 .  Nume  r ic  a  1  Prob 1 ems 

The  covariance  matrix  P  at  time  t  for  all  of  the 
elemental  filters  proved  to  be  numerically  ill-conditioned 
(a  function  of  there  being  a  accurate  data  on  some  of  the 
states  from  the  measurements).  Therefore  it  was  impossible 
to  compute  the  ambiguity  function  as  described  by  Equation 
(11-33),  as  that  equation  requires  the  determinant  of 
P(t£+).  Since  the  P(t^+)  matrix  was  used  to  compute  P(t^  ), 
which  is  used  in  the  computation  of  the  A  matrix  used  in 
Equations  (II-7),  (11-25),  (11-30)  and  (11-33),  more 

problems  resulted,  as  the  computed  determinants  of  all  of 


the  A  matrices  proved  to  be  negative.  Tunning  of  the 
filters  could  result  in  improving  the  numerical  condition  of 
PCt^-).  and  could  prevent  the  A  matrices  determinants  from 
being  negative.  However;  this  would  envolve  the  use  of  less 
realistic  noise  strengths  for  Q  and  R,  this  could  be 
worthwile  as  an  academic  exercise  but  time  constraints 
prevented  that  from  being  done  for  this  thesis. 

In  order  to  overcome  these  numerical  difficulties,  the 
expressions  for  the  probability  weighting  factors  (Equation 
(11-25)),  and  the  ambiguity  function  (Equation  (11-33)  were 
approximated  to  remove  the  necessity  to  use  the  determinants 
of  the  P(t^*)  and  the  A  matrices.  In  Equation  (11-25) 
the  density  function  is  approximated  as: 

fj(z(ti))  -  exp[-(l/2)£jT(ti)Aj*1(ti)rj(t.)] 

Note  this  is  no  longer  a  true  density  function  because  the 
scale  factor  is  now  incorrect;  but  because  of  the 
denominator  in  Equation  (11-25)  the  probability  weightings 
are  still  scaled  correctly.  This  change  is  also  reflected 
in  Equation  (11-30)  describing  how  the  probability  weighting 
for  a  new  filter  in  the  bank  can  be  computed  based  on  its 
expected  correctness.  This  approximation  is  reasonable  as 
long  as  the  A  matrices  of  the  the  elemental  filters  are 
close  enough  to  each  other  so  that  their  determinants,  in 
the  absence  of  numerical  problems,  would  be  expected  to  be 
of  approximately  equal  magnitudes.  The  ambiguity  function 
described  in  Equation  (11-33)  was  approximated  by  removing 


the  terms  containing  the  determinants  of  P(ti  +  )  and  A ;  it 
then  became: 

Ai(A.At)  “  m/2  ln(2n)  -  n/2  ln(2n) 

-  1/2  tr{A_1(ti;A)  [H(ti)Pe(ti";At*a)HT(ti)+R(ti)]} 

-  1/2  tr{P“1(ti+;A)Pe(ti+;at,A)} 

This  approzc imat ion  is  reasonable  as  the  determinants  of 
P(t£+)»  and  A  (in  the  absense  of  the  numerical  difficulties) 
can  be  expected  to  have  minimal  impact  on  the  ambiguity 
function  as  the  primary  sensitivity  of  the  ambiguity 
function  is  in  the  quadratic  terms  which  are  being 
preserved . 

11.3  Amb iauitv  Fftaslifla  Analysis 

Analysis  of  the  ambiguity  function  was  based  on  three- 
dimensional  view  plots  which  show  the  magnitude  of  the 
ambiguity  function  as  a  continuous  surface  over  the 
two  dimensional  parameter  space,  for  specific  true  parameter 
values.  Recall  that  the  two  uncertain  parameters  are  the 
mass  density  of  the  spacecraft  arms  and  the  modulus  of 
elasticity  of  the  arms.  Three  typical  plots  are  included  as 
Figures  V-l,  V-2,  and  V-3,  based  on  true  parameter  points 
(5,5),  (3,7),  and  (7,3)  respectively,  the  arrow  on  the  plot 
indicates  the  true  parameter  point.  The  numbers  for  the 
uncertain  parameters  are  not  the  true  values  of  the 
parameters  at  those  points,  but  indices  indicating  which  of 
the  ten  discrete  values  for  each  uncertain  parameter  is 
used.  As  can  be  seen  the  ambiguity  functions  are  fairly 


flat,  which  indicates  possible  sensitivity  problems,  and 
there  are  several  ripples,  which  indicates  possible  problems 
identifying  the  correct  parameter  values.  The  large  peak  at 
the  (0,0)  point  should  be  ignored  as  it  is  an  artificial 
result  of  the  curve  fitting  routine  used  to  generate  the 
contours  and  of  the  fact  that  there  is  no  data  point  at 
(0,0).  In  addition,  the  plots  for  different  true  parameter 
points  are  very  similar.  This  indicates  that  the  moving 
bank  multiple  model  adaptive  estimator  will  have 
difficulties  identifying  the  correct  parameter  point,  this 
is  especially  true  when  the  true  parameter  point  is  not 
located  on  at  least  a  local  peak,  as  is  the  case  in  Figure 
V-3  for  true  parameter  point  (7,3).  However,  even  if  the 
true  parameter  point  cannot  be  identified,  the  state 
estimates  may  still  be  accurate,  if  the  different  models  are 
all  doing  a  good  job  of  estimating  the  states  even  though 
based  on  the  incorrect  system  model.  The  true  parameter 
point  not  being  located  on  a  local  peak  in  the  ambiguity 
surface  as  occurs  in  Figure  V-3  was  not  anticipated,  nor 
believed  possible,  and  is  likely  due  to  the  approximations 
used  in  the  calculation  of  the  ambiguity  function.  However, 
a  rigorous  proof  that  the  true  parameter  point  must  lie  on 
at  least  a  local  peak  is  not  easily  done  and  is  beyond  the 
scope  of  this  thesis. 

The  ripples  and  multiple  local  peaking  argue  for  the 
use  of  a  contraction  algorithm  to  aid  initial  acquisition  of 
the  true  parameter  point.  Without  a  contraction  algorithm. 


the  moving  bank  mast  start  at  its  smallest  discretization 


and  move  to  center  itself  on  the  local  peak  with  the 
greatest  height.  If  a  local  peak  not  due  to  the  true 
parameter  point  is  within  the  area  of  the  bank,  or  if  while 
the  bank  moves  it  reaches  an  erroneous  peak  first,  the  bank 
may  remain  there  (noise  may  make  it  move  off).  The  goal  of 
the  contraction  algorithm  is  not  only  to  speed  the  initial 
movement  of  the  bank,  but  by  covering  more  area,  the  bank 
may  be  able  to  contract  to  the  neighborhood  of  the  correct 
parameter  point  and  ignore  local  peaks  not  doe  to  the  true 
parameter  point. 

Conversely,  the  flatness  of  the  surfaces  argues  against 
residual  monitoring  using  the  likelihood  quotient  being  able 
to  signal  an  expansion  of  the  bank  by  identifying  when  a 
jump  change  in  the  true  parameter  point  has  taken  place.  On 
a  flat  ambiguity  function  surface  it  is  likely  that  changes 
in  the  residuals  due  to  the  jump  change  will  be  masked  by 
the  normal  changes  in  residuals  seen  as  the  effects  of  the 
dynamics  and  measurement  noises.  For  residual  monitoring 
using  the  likelihood  quotient  to  work,  the  jump  change  in 
the  true  parameter  point  has  to  cause  a  change  in  the 
residuals  greater  than  the  flucuations  due  to  noise.  If  not 
a  jump  event  detection  threshold  level  can  not  be  set, 
because  the  noise  changes  would  signal  a  jump  change. 

The  flatness  of  the  ambiguity  surfaces  also  indicates  that 
adaptive  estimation  and  control,  of  any  kind,  is  probably 
not  worthwhile  for  this  system.  In  fact,  a  filter  based  on 
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a  nominal  parameter  point  may  provide  adequate  state 
estimates  and  an  acceptable  controller  gain  matrix.  Of 


course  the  definition  of  nadequate  state  estimatesn,  and 
nacceptable  controller  gain  matrixn  are  based  on  the  level 
of  accuracy  required  by  the  individual  implementation.  The 
flatness  of  the  surfaces  was  not  anticipated,  considering 
the  large  changes  in  the  parameters  used  to  create  the 
parameter  space,  and  the  core spond ing ly  large  changes  in  the 
eigenvalues  of  the  bending  modes  that  were  observed  (see 
Appendix  A  for  a  listing  of  the  eigenvalues).  The 
approximations  used  in  computing  the  ambiguity  function  may 
be  one  reason  that  the  surfaces  are  so  flat.  The  wealth  of 
measurement  data  (five  position  measurements  and  five 
acceleration  measurement)  available  for  this  model  also 
worked  to  minimize  the  impact  of  the  changes  in  the  model 
used  to  set  up  the  paramater  space  upon  the  precision  of 
state  estimation  or  control.  This  was  seen  by  examining  the 
Kalman  filter  gain  matrices  in  conjunction  with  the 
measurement  matrices  and  observing  that  a  heavy  weight  was 
put  on  the  measurements  as  opposed  to  the  estimated  states 
before  the  measurements.  The  inaccurate  estimated  states 
are  thus  essentially  ignored. 

The  academic  solution  to  the  difficulties  encountered 
at  this  point  might  be  to  reduce  the  number  or  precision  of 
the  measurements  to  obtain  a  problem  for  which  adaptive 
estimation  and  control  would  be  more  useful.  For  this 
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problem  it  would  be  more  useful  to  reduce  the  amount  of 
adaptivity  by  only  using  one  uncertain  parameter,  or  by 
taking  the  adaption  off  line  entirely  and  only  using  a 
robustified  f i 1 t e r / con t r o 1 1 e r.  However;  despite  the 
potential  difficulties  forcast  by  the  ambiguity  function 
analysis,  the  addaptive  estimator  was  applied  to  this 
problem.  This  was  done  to  confirm  the  usefulness  of  the 
ambiguity  function  as  a  predictor  of  estimator  performance, 
and  to  evaluate  how  large  an  effect  the  difficulties  had  on 
estimator  perfomance.  The  following  sections  describe  the 
algorithm  characteristics  for  this  application. 

V . 4 .  Bank  Movement 

The  first  simulations  started  the  moving  bank  at  it 
finest  discretization  level,  and  located  the  bank  center  at 
(5,5),  chosen  as  the  approximate  center  of  the  parameter 
space,  and  allowed  it  to  move  using  the  probability 
monitoring  movement  algorithm  developed  by  Maybeck  and  Hentz 
[26,27]  described  in  Chapter  II,  where  the  bank  is  centered 
over  the  filter  having  the  largest  probability  weighting. 

The  points  (2,2),  (3,7),  (5,5),  (8,8),  and  (7,3)  where 
chosen  as  the  true  parameter  points  to  test  the  bank's 
ability  to  move  in  different  directions.  Before  making 
multiple  Monte  Carlo  runs  to  test  the  probability 
initialization  algorithms,  single  runs  were  made  to  set  the 
probability  move  threshold,  and  to  determine  a  strength  for 
the  dither  signal  be  used  to  excite  the  system  states  and 


allow  parameter  point  acquisition  in  the  absence  of  any 
feedback  control.  A  dither  signal  of  200  amplitude,  divided 
equally  between  the  three  inputs,  at  20  Hz  was  applied  to 
all  the  control  inputs.  The  dither  signal  was  determined  by 
trial  and  error  as  a  determination  of  an  optimal  signal  (see 
[  5  : 68-15 1 , 12 : 223-260 ]  for  a  discussion  of  optimal  inputs) 
would  be  beyond  the  scope  of  this  thesis. 

Maybeck  and  Hentz  found  that  a  probability  move 
threshold  of  0.05  enabled  the  bank  to  move  anytime  the  bank 
was  not  centered  on  the  filter  having  the  highest 
probability  weighting  [26:52,27:19].  This  will  be  true  for 
any  threshold  below  1/9  =  0.111,  which  is  the  lowest 
probability  weighting  the  center  filter  can  have,  and  still 
be  weighted  equal  to,  or  heavier  than  the  other  filters.  In 
this  application  however,  it  was  found  that  a  threshold 
slightly  larger  than  0.111  speeded  acquisition  of  the  true 
parameter  point  by  eliminating  movement  (possibly  in  the 
wrong  direction)  due  to  noise  when  the  filter  probability 
weightings  were  all  close  in  value.  A  threshold  of  0.15 
allowed  the  bank  to  center  itself  on  the  (3,7)  true 
parameter  point  after  0.98  seconds  vs.  1.54  seconds  with  a 
threshold  of  0.05.  Increasing  the  threshold  further  to  0.2 
however  increased  the  acquisition  time  to  1.28  seconds.  A 
threshold  level  of  0.15  was  therefore  used  for  all 
subsequent  runs.  For  the  single  runs,  the  p s e ud o- r andom 
number  generator  used  to  generate  the  noise  signals  was 
always  seeded  with  the  same  starting  number  to  ensure  that 


the  runs  were  directly  comparable 


Monte  Carlo  studies  of  the  bank's  ability  to  move  to 
each  of  the  five  chosen  true  parameter  points,  using  ten 
runs  apiece,  were  made  for  both  the  simple  method  of 
initializing  the  probability  weightings  of  the  new  filters 
brought  into  the  bank  by  weighting  each  new  filter  equally, 
and  for  the  more  complicated  method  of  estimating  a 
probability  weighting  for  each  new  filter  based  on  Equation 
(11-30).  The  results  can  be  seen  in  Appendix  B  in  Figures 
B-l  through  B-10. 

These  figures  are  plots  of  the  statistics  (mean  +  one 
standard  deviation)  of  the  error  in  the  parameter  estimates. 
EA(1)  is  the  error  in  the  estimate  of  the  parameter 
reflecting  the  mass  density  of  the  arms,  and  EA(2)  is  the 
error  in  the  estimate  of  the  parameter  reflecting  the 
modulus  of  elasticity  of  the  arms.  Once  the  matrices 
describing  the  system  at  a  particular  parameter  point  have 
been  calculated,  the  true  values  of  the  two  uncertain 
parameter  were  no  longer  used.  It  is  far  more  convienent  to 
index  the  parameter  values  to  the  whole  numbers  one  through 
ten  reflecting  their  position  in  the  parameter  space.  Thus, 
the  (1,1)  parameter  point  is  where  the  uncertain  parameters 
take  on  the  discrete  value  20%  less  than  their  nominal 
values,  and  the  point  (6,6)  is  where  they  are  equal  to  their 
nominal  values.  An  increase  of  one  in  the  parameter  index 
corresponds  therefore,  to  an  increase  of  4%  of  the  nominal 
value.  The  errors  in  the  parameter  estimates  are  calculated 


in  terms  of  the  parameter  indices;  if  the  true  parameter 
point  is  at  (4,5)  and  the  location  of  the  estimated 
parameter  vector  is  (3. 5, 6.0),  the  error  in  the  parameter 
estimate  plots  will  show  a  value  of  0.5  for  the  mass  density 
parameter,  and  a  value  of  -1.0  for  modulus  of  elasticity 
parameter.  It  is  important  to  note  that  because  the 
estimated  parameter  vector  is  a  weighted  average  of  the 
indices  of  the  discrete  parameter  points  whose  associated 
filters  are  currently  implemented  in  the  moving  bank,  the 
location  of  the  estimated  parameter  vector  is  not  likely  to 
be  exactly  equal  to  one  of  the  discrete  parameter  points. 

The  first  set  of  plots  shows  the  error  in  the  parameter 

estimate  («a)  when  the  probability  weightings  of  the  new 
filters  in  the  bank  were  calculated  using  Equation  (11-30) 
(Figures  B-l  through  B-5).  The  second  set  of  plots  shows  s.a 
when  the  new  filters  are  initialized  by  equally  distributing 
the  unused  probability  (Figures  B-6  through  B-10).  A 
comparison  shows  that  there  is  no  appreciable  difference  in 
the  results  of  the  two  methods.  However,  it  is  a 
considerably  larger  computational  burden  (by  several  orders 
of  magnitude,  depending  on  the  architecture  of  the 
processor)  to  compute  the  probability  weightings  using 
Equation  (11-30). 

The  performance  predicted  by  the  ambiguity  function 

analysis  can  also  be  seen  in  the  plots  of  £a  error 
statistics  for  the  case  of  only  moving  the  bank  without 
changing  its  size,  using  either  of  the  filter  initialization 


algorithms.  In  all  cases  the  estimator  does  a  better  job  of 
estimating  the  second  parameter  than  the  first:  it  is  always 
identified  within  one  standard  deviation  by  the  time  one 
second  has  elapsed.  However,  only  the  true  parameter  points 
at  (5,5)  and  (3,7)  were  correctly  identified,  for  both 
parameters,  within  one  standard  deviation.  These  statistics 
are  matched  by  center  of  the  bank  statistics:  in  all  cases 
the  bank  center  had  moved  to  the  correct  second  parameter  by 
0.8  seconds. 

Examination  of  the  ambiguity  function  plots  shows 
ridges  running  diagonally  across  the  surfaces.  The  bank 
will  not  move  down  into  the  valleys  between  the  ridges, 
seeking  instead  the  peaking  regions  of  the  abiguity  function 
surface.  Therefore,  when  the  true  parameter  point  is  not  at 
a  local  peak  accessable  from  the  ridge  on  which  the  starting 
point  (5.5)  is  located,  the  estimator  is  forced  to  make  a 
choice  as  to  w  h  i  c  h  p  a  r  a  m  e  t  e  r  it  will  do  a  good  job 
estimating.  The  reason  it  is  the  second  parameter  that  is 
favored  can  be  seen  by  examining  Tables  A  - 1  and  A  -  2  in 
Appendix  A.  These  tables  show  the  eigenvalues  of  the  second 
and  forth  bending  modes  at  the  various  parameter  points.  It 
can  be  seen  that  a  change  in  the  second  parameter  (E;  the 
modulus  of  elasticity)  causes  a  larger  change  in  the 
eigenvalues  than  does  a  change  in  the  first  parameter  (p: 
the  mass  density).  Thus  the  filters  on  the  ridge  with  the 
more  correct  value  of  E  recieve  a  higher  probability 
weighting  than  the  filters  with  the  more  correct  value  of 
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density,  and  the  bank  moves  toward  the  correct  value  for  the 
second  parameter. 

The  overall  flatness  of  the  ambiguity  surfaces  was 
reflected  in  the  plots  of  the  statistecs  of  the  errors  in 
the  state  estimates  <*x).  These  plots  showed  ex  to  be 

essentially  independent  of  e«.  Figure  B-ll  shows  &x 
statistcs  for  the  case  where  the  true  parameter  point  was  at 

(3,7),  and  Figure  B-12  depicts  &x  for  th«  case  where  the 
true  parameter  point  was  at  (8,8).  The  states  1-6  reflect 
the  position  of  the  rigid  body,  the  second  bending  mode,  and 
the  fourth  bending  mode  respectively,  and  then  the 
velocities  in  the  same  order.  Both  sets  of  plots  are 
identical,  as  were  the  plots  no  matter  where  the  true 
parameter  point  was  set.  The  point  is  especially 
illustrated  by  Figures  B-13,  and  B-14,  which  show  ex 
characteristics  for  just  a  single  filter  based  on  assumed 
parameter  values  at  (5,3).  In  the  first  case  (Figure  B-13) 
the  true  parameter  point  was  at  (3,7),  and  in  the  second 
case  the  true  parameter  point  was  actually  at  (5,5):  the 
filter  with  parameters  at  (5,5)  does  an  equally  good  job  of 
estimating  the  states  in  both  cases. 

V . 5 .  Bank  Contraction 

Begining  the  run  with  an  acquisition  phase  using  a 
coarser  bank  discretization,  then  contracting  the  bank  to 
the  finest  discretization  in  two  steps  was  tested  to  see  if 
it  speeded  parameter  acquisition.  Parameter  error 


covariance  monitoring  was  used  to  detect  the  acquisition  of 


the  parameters,  and  signal  the  change  in  discretization. 

The  scalar  measures  associated  with  the  parameter  estimate 
covariance  matrix,  selected  to  signal  contraction,  were  the 
parameter  variances  (the  diagonal  elements  of  the  covariance 
matrix).  The  variances  of  both  parameters  were  compared  to 
a  contraction  threshold  and  the  result  was  used  to  signal 
contract  ion. 

Single  runs  were  used  to  set  the  contraction 
thresholds.  Two  were  required,  one  for  each  level  of 
contraction.  The  first  coarsest  discretization  level  places 
the  filtesr  four  parameter  points  apart.  The  second 
intermediate  discretization  level  places  the  filters  two 
parameter  points  apart.  And  the  final  finest 
discretizat ionlevel  places  the  filters  one  parameter  point 
apart.  An  additional  consideration  was  whether  to  require 
both  or  only  one  of  the  variances  to  be  below  the  threshold 
in  order  to  signal  the  contraction.  It  was  found  that 
requiring  both  variances  to  be  below  the  threshold  gave  a 
more  accurate  determination  of  when  contraction  was 
warranted;  however,  this  required  tradeoffs  in  the 
determination  of  the  threshold,  as  the  same  threshold  was 
not  necessarily  most  appropriate  for  both  variances.  The 
threshold  values  were  then  chosen  equal  to  the  larger  of  the 
two  variances  at  the  time  when  the  true  parameter  point 
would  lie  within  the  area  of  the  bank  after  it  contracted 
about  the  estimated  parameter  point.  This  required 
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artificial  knowledge  of  the  true  parameter  point,  and  the 
printout  of  the  time  period  by  time  period  statistics  on  the 
bank  location  and  estimates.  A  threshold  was  thns  selected 
for  each  of  the  four  test  points  (the  same  test  points  as 
for  the  move  analysis  with  the  exception  of  ( 5,5 ))  and 
smallest  of  the  five  then  selected  as  the  best  overall. 

This  threshold  level  setting  mechanism  worked  well  for 
signaling  contraction  from  coarsest  discretization  level  to 
the  intermediate  discretization  level.  Bnt  it  conli  not  be 
used  in  all  cases  to  select  a  threshold  for  signaling  the 
second  contraction  to  the  finest  discretization  level.  When 
the  true  parameter  point  was  at  (8,8),  the  parameter 
estimate  never  moved  closer  to  the  true  parameter  point 
after  the  first  contraction.  Therefore  a  contraction 
threshold  based  only  on  the  other  three  points  was  used. 

The  contraction  thresholds  chosen  were  6.0  for  the  four-to- 
two  contraction,  and  2.6  for  the  two-to-one  contraction. 

The  contraction  thresholds  proved  to  be  sensitive  to 
the  probability  weight  lower  bound  used  to  keep  the  bank 
from  locking  on  to  a  single  filter.  When  the  probability 
weight  lower  bound  is  set  too  high,  it  prevents  the  the 
variances  from  decreasing  below  the  thresholds  in  all  cases. 
Therfore.  the  contraction  thresholds  need  to  be  reset  if  the 
probability  lower  bounds  change.  Lower  bounding  the 
probability  weights  also  causes  a  larger  error  in  the 
parameter  estimate,  and  possibly  the  state  estimate 
(although  that  was  not  true  in  this  implementation),  as  it 
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requires  an  artificially  high  probability  weight  be  assigned 


to  a  poor  estimate  in  steady  state.  The  purpose  in  lower 


bounding  the  probability  weights  is  to  allow  the  bank  to 


react  to  a  change  in  the  true  parameter  point  after  the 


initial  true  parameter  point  has  been  identified.  There  is 


thus  a  tradeoff  to  be  made  between  keeping  the  probability 


lower  bound  high  enough  to  enable  quick  reaction  to  changes 


in  the  true  parameter  point,  and  low  enough  not  to  interfere 


with  the  contraction  algorithm  nor  cause  too  great  an  error 


in  the  parameter  and  state  estimates.  A  lower  bound  of  0.01 


worked  well  in  this  implementation;  when  it  was  raised  to 


0.05,  the  bank  would  not  contract  from  the  coarsest 


discretization. 


The  results  of  using  the  contraction  algorithm  can  be 


seen  in  Figures  B-15  through  B-18,  which  show  the  plots  of 


£a  statistics  vs.  time  for  the  true  parameter  points  (2,2), 


(3,7),  (7,3),  and  (8,8)  respectively.  The  results  can  be 


compared  to  the  plots  of  £.a  characteristics  vs.  time  for  the 


movement-only  options  depicted  in  Figures  B-l  through  B-10. 


It  can  be  seen  that  the  use  of  the  contraction  aquisition 


phase  is  successful  in  correctly  identifying  parameter  two 


(E)  more  quickly;  however  only  small  improvements  in 


reducing  the  error  in  parameter  one  (p)  are  achieved  and 


only  temporarily.  The  expense  is  the  larger  standard 


deviations  which  result  from  the  coarser  discretizations. 


In  terms  of  the  shape  of  the  ambiguity  function;  the 


acquisition  cycle  is  successful  because  the  flatness  of  the 
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surface,  and  the  ripples  and  local  peaks  slow  the  movement 
of  the  bank. 

There  is  very  little  added  (and  perhaps  less) 
computation  due  to  the  contraction  algorithm.  Either  one  or 
three  additional  comparisons  are  required  in  the  logic  (one 
to  determine  if  the  bank  is  at  a  coarse  discretization;  if 
so,  each  of  the  two  parameter  error  variances  is  compared  to 
the  appropriate  threshold),  but  the  operations  required  for 
the  actual  contraction  are  very  similar  to  those  required 
for  a  move  and,  as  all  of  the  probability  weights  are  reset 
to  0.111  upon  contraction,  less  initialization  computation 
is  required  than  is  required  following  a  bank  move.  In 
addition  some  computation  is  saved,  since  with  contraction 
fewer  moves  are  required  to  get  the  bank  to  its  final 
location.  Thus  computation  is  saved  in  the  short  run,  but 
as  the  run  continues,  in  the  absense  of  a  re-expansion  to 
the  coarse  discretization,  the  total  computation  will 
eventually  be  slightly  greater.  If  re-expansion  occurs 
following  the  signaling  of  a  jump  change  in  the  true 
parameter  values  there  is  again  a  savings  during  the 
r e acqu i s  i  t ion  of  the  true  parameter  values.  The  real  issue 
here  is  the  amount  of  computation  per  sample  period;  and  the 
contraction  algorithm  is  essentially  equal  to  the  non¬ 
contraction  algorithm  in  that  regard. 


V .  6 . 
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It  was  hoped  that  residual  monitoring  would  be  able  to 
detect  jump  changes  in  the  true  parameter  point.  Unfortu¬ 
nately,  when  an  attempt  to  set  a  detection  threshold  was 
made,  the  likelihood  function  proved  much  more  sensitive  to 
the  noise  inputs  than  to  the  difference  in  the  residuals  due 
to  the  change  in  the  true  parameter  point.  Vhen  the  attempt 
to  set  a  threshold  was  made,  on  one  run  the  likelihood 
function  actually  decreased  at  the  time  of  the  jump  change 
due  to  a  smaller  noise  input  that  sample  period.  The 
failure  of  the  jump  change  detection  by  residual  monitoring 
can  be  explained  in  terms  of  the  flatness  of  the  ambiguity 
function  as  noted  earlier.  Therefore,  no  evaluation  of  the 
effectiveness  of  a  reacquisition  period,  using  the  bank 
expansion  and  contraction  algorithms,  following  a  jump 
change  in  the  true  parameter  point  could  be  made. 

A  plot  of  ea  statistics  vs.  time  for  a  case  where  the 
bank  was  started  at  its  largest  discretization,  allowed  to 
contract  and  move  to  a  true  parameter  point  at  (3,7),  and 
react  to  a  jump  change  in  the  parameter  point  to  (5,5)  at 
1.0  second,  reveals  that  the  move  algorithm  does  continue  to 
track  the  true  parameter  point  via  bank  movement  without  a 
change  in  the  bank  discretization  level  (Figure  B-19).  As 
before  parameter  two  is  consistently  estimated  better,  both 
before  and  after  the  jump  change.  For  the  same  situation. 
Figure  B-20  reveals  that,  as  before,  the  state  estimates  for 
this  implementation  do  not  suffer. 
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V .  7  Coa.i£0llex  s 

The  evaluation  of  the  effectivness  of  the  various 
adaptive  control  algorithms  vas  also  abandoned  due  to  the 
accuracy  of  the  state  estimates,  independent  of  the  accuracy 
of  the  parameter  estimate.  In  addition,  inspection  of 

sample  matrices  reveals  little  change  due  to  p,  and 
almost  none  due  to  E;  so  no  large  errors  in  the  control 
input  can  be  expected  as  a  result  of  the  use  of  the  wrong 
matrix . 

V . 8  Summary 

The  results  of  the  simulations  were  presented  in  this 
Chapter.  Initial  problems  due  to  numerical  difficulties 
were  overcome  through  the  use  of  judicious  aproximat  ions. 
Analysis  of  the  ambiguity  function  revealed  several  insights 
into  the  expected  performance  of  the  multiple  model  adaptive 
estimator,  and  the  probable  lack  of  a  need  for  an  adaptive 
controller.  These  insights  were  borne  out  by  the  actual 
simulation  results.  The  two  movement  algorithms  worked 
equally  well  but  were  unable  to  move  across  the  valleys 
depicted  in  the  ambiguity  function  plot.  The  estimator 
proved  more  sensitive  to  the  second  uncertain  parameter  (E  = 
the  modulus  of  elasticity  of  the  arms)  than  to  the  first  (p 
-  the  mass  density  of  the  arms).  The  contraction  algorithm 
gave  some  increase  in  performance,  but  essentially  only  for 
the  second  uncertain  parameter.  Importantly,  none  of  the 
movement  or  contraction  logics  affected  the  overall  accuracy 


of  the  state  estimates  in  any  discernible  manner.  Residual 
monitoring  proved  ineffectual  in  detecting  jump  changes  in 
true  parameter  points,  making  an  evaluation  of  a 
reacquisition  phase,  using  the  bank  contraction  and 
expansion  algorithms,  impossible.  The  evaluation  of  the 
adaptive  controllers  was  also  abandoned  due  to  the  accuracy 
of  the  state  estimates.  Conclusions  regarding  these  results 
and  recommendations  for  areas  of  further  study  are  presented 
in  Chapter  VI. 
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Ambiguity  function  analysis  was  able  to  predict/explain 
most  of  the  problems  encountered  in  the  simulations  of  the 
multiple  model  adaptive  estimator/controller.  It  proved  to 
be  a  valuable  tool  that  should  be  used  early  in  the  design 
of  the  adaptive  estimator  in  order  to  aid  in  the  selection 
of  appropriate  parameters  for  estimation,  and  of  workable 
move/contract/expansion  stategies.  This  is  especially  true, 
as  the  results  of  the  simulations  showed  that  each  of  the 
strategies  investigated  have  shortcomings  which  need  to  be 
overcome.  The  bank  may  not  be  able  to  reach  all  points  by 
movement  alone,  or  at  all.  This  occurs  when  the  points  are 
located  in  a  valley  in  the  ambiguity  function,  or  if  the 
bank  encounters  a  local  peak  while  moving  towards  the 
correct  parameter  values.  Square  contraction  may  require 
compromise  in  the  setting  of  the  contraction  thresholds, 
which  will  hinder  its  performance.  This  occurs  when  the 
bank  does  a  better  job  of  estimating  one  parameter  than 
another.  Separate  contraction  thresholds  which  enable  the 
bank  to  contract  separately  in  each  direction,  when  it  is 
most  appropriate  to  contract  in  that  direction,  should 
improve  performance.  And  residual  monitoring  via  the 
likelihood  quotient  may  not  be  able  to  detect  jump  changes, 
depending  on  the  relative  effect  on  the  residuals  of  noise 


vs.  system  changes. 


VI. 2.1  Ambiguity  funct ion  Analysis.  The  first 
question  that  early  ambiguity  function  analysis  can  answer 
is  whether  adaptive  estimation/control  is  needed  at  all. 

This  is  not  always  known  in  the  initial  stages  of  system 
design,  as  was  the  case  for  the  model  used  for  this 
investigation.  Once  it  is  determined  that  adaptive 
e s t imat ion/ control  is  required,  the  ambiguity  function 
analysis  can  determine  which  of  the  uncertain  parameters 
need  to  be  estimated  and  which  maybe  ignored.  It  maybe 
possible  to  combine  multiple  uncertain  parameters  into  a 
single  artificial  parameter,  or  into  a  reduced  number  of 
artificial  parameters,  to  be  used  by  the  adaptive 
estimator/controller.  This  enables  the  estimator  to 
concentrate  on  estimating  the  uncertain  system  parameters 
(either  physical  quantities,  or  deriv.d  quantities)  that 
have  the  greatest  impact  on  system  performance. 

Once  the  parameter  space  is  chosen,  further  analysis  of 
the  ambiguity  function  may  be  useful  in  chosing  the 
discretizations  used  as  a  basis  for  elemental  filter  design. 
The  discretization  could  be  finer  in  regions  of  greater 
parameter  variation  sensitivity,  and  coarser  in  regions  of 
lesser  parameter  variation  sensitivity,  instead  of  being 
spread  evenly  over  the  entire  parameter  space.  Ambiguity 
function  analysis  may  then  be  helpful  in  chosing  the  initial 
conditions  for  the  moving  bank:  bank  center  and 
discretization  level.  As  seen,  from  certain  starting 
positions,  the  bank  may  not  be  able  to  move  to  the  true 


parameter  point.  Judicious  choice  of  the  starting  position 
may  alleviate  this  problem  for  a  large  class  of  true 
parameter  points.  Then  thought  can  be  given  to  the  bank 
contraction  and  ezpandsion  strategies  to  be  implemented.  If 
multiple  local  peaks  are  present,  contraction  strategies  may 
be  especially  valuable,  as  contraction  about  the  estimated 
parameter  values  could  prevent  the  bank  from  moving  to,  and 
remaining  on,  a  local  peak  positioned  between  the  bank 
starting  position  and  the  true  parameter  point. 

VI. 2. 2.  Bank  Movement.  Despite  the  overall  lack  of  a 
requirement  for  adaptive  estimation  in  this  application. 


because  all  of  the  elemental  filters  did  an  acceptable  job 
of  state  estimation,  the  probability  monitoring  method  of 
moving  the  bank  was  able  to  do  a  surprisingly  good  job. 

Even  when  the  bank  was  not  able  to  center  itself  over  the 
true  parameter  point,  it  was  able  to  move  in  the  right 
general  direction.  Two  refinements  to  the  algorithm  were 
discovered.  Calculation  of  the  probability  weightings  for 
initialization  of  new  filters  in  the  bank  based  on  the 
distance  from  the  estimated  true  parameter  point  is  not 
worthwhile.  Evenly  deviding  the  total  probability  weighting 
taken  from  the  filters  removed  from  the  bank,  amongst  those 
added  to  the  bank  gave  equally  good  performance  with  less 
computation.  And,  it  may  be  worthwhile  to  raise  the  bank 
movement  threshold  above  0.111,  in  order  to  prevent  noise 
jitter  movement  of  the  bank.  If  the  threshold  is  set  below 
0.111,  it  is  unnecessary  and  can  be  eliminated,  simply  check 
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to  insure  the  bank  is  centered  on  the  filter  with  the 
highest  probability  weighting. 

VI. 2. 3.  Contraction.  Contraction  of  the  bank  from  an 
initial  coarse  discretization  to  a  fine  discretization, 
proved  to  be  nsefnl  in  improving  the  speed  of  parameter 
acquisition.  Little  additional  computation  is  involved  in  the 
use  of  this  strategy;  however,  there  is  a  price  to  pay  in 
larger  standard  deviations  of  the  parameter  estimate  error, 
as  noise  causes  a  greater  change  in  the  parameter  estimate 
when  the  bank  is  in  the  coarser  discretization. 

The  best  contraction  thresholds  for  the  two  uncertain 
parameters  were  not  identical  and  a  compromise  threshold  had 
to  be  chosen.  This  meant  that  contraction  did  not  always 
move  the  bank  center  in  the  right  direction  relative  to  the 
true  value  of  the  first  uncertain  parameter,  wheras  the 
bank  center  always  move  towards  the  true  value  of  the  second 
uncertain  parameter.  A  two-step  contraction  was  used,  but 
an  advantage  of  the  two  step  contraction  vs.  a  one-step 
contraction  is  unproven. 

VI. 2. 4.  Res idual  Monitoring.  Monitoring  the  residuals 
through  the  use  of  a  likelihood  quotient  in  order  to  detect 
a  jump  change  in  the  true  parameter  vector  value  did  not 
prove  useful  because  of  the  lack  of  sensitivity  in  either 
state  estimate  precision,  or  residual  characteristics,  to 
assumed  paremeter  value.  The  magnitude  of  the  likelihood 
quotient  varied  more  as  a  result  of  the  system  noises  than 
as  a  result  of  the  changes  in  the  uncertain  parameters, 
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making  it  impossible  to  select  a  threshold  valne  at  which  to 
trigger  an  expansion.  It  must  be  remembered  however,  that 
the  residuals  are  the  only  real  source  of  information  on  the 
trne  system  and  mnst  be  the  basis  for  all  move / cont r ac t / 
expand  strategies. 

VI. 3.  Areas  of  Far  the r  Stndv 

Further  study  in  all  areas  of  moving  bank  multiple  model 
adaptive  control  is  still  needed.  A  specific  area  where 
further  study  is  needed  was  revealed  by  the  importance 
ambiguity  function  analysis  was  found  to  have.  If  the 
ambiguity  function  could  be  shaped,  by  changes  in  the 
uncertain  parameters  chosen  to  be  estimated,  or  by  changes  in 
the  uncertain  parameter  discretization,  to  reflect  desirable 
characteristics  for  the  move/ con trac t/ expand  logics, 
e s t i ma to r / con t r o 1 1 er  performance  may  be  enhanced.  The 
ambiguity  function  could  be  shaped  through  selection  of 
artificial  uncertain  parameters  and  choice  of  discretized 
parameter  points. 

An  evaluation  of  the  adv an t age s / d  i  s advan t age s  of 
rectangular  contraction  vs.  square  contraction  to  eliminate 
the  undesirable  compromises  necessary  in  the  selection  of  a 
contraction  threshold  would  be  worthwhile.  Rectangular 
contraction  could  allow  the  contraction  in  each  direction 
when  it  is  most  appropriate,  making  bank  size  in  one 
direction  independent  of  bank  size  in  the  other. 

Rectangular  contraction  could  prevent  the  movement  of  the 


bank  center  away  from  the  true  value  of  one  of  the  uncertain 
parameters  because  of  the  current  inacurate  estimate  of  that 
parameter.  These  posible  advantages  would  have  to  be 
weighed  against  the  required  increase  in  computational 
loading.  Evaluation  of  two-step  vs.  one-step  contraction 
under  different  conditions  might  reveal  when  one  is  more 
valuable  than  the  other.  Other  methods  of  triggering  the 
contractions  could  also  be  workable.  It  may  be  possible  to 
contract  away  from  filters  with  very  small  probability 
weightings. 

Methods  of  detecting  jump  changes  in  the  uncertain 
parameter  values  still  need  further  investigation.  Residual 
monitoring  using  the  likelihood  function  may  be  adequate  in 
aplications.  unlike  the  application  investigated  here,  where 
wrongly  assumed  parameter  values  do  affect  state  estimation. 
It  may  be  possible  to  make  modifications  to  the  likelihood 
quotient  method  to  prevent  noise  from  masking  the  results  of 
uncertain  parameter  jump  changes  by  adding  additional 
conditions.  An  evaluation  of  the  changes  in  the  probability 
weightings  relative  to  each  other  may  also  be  capable  of 
indicating  jump  changes.  If  all  of  the  probabilities 
converge  on  a  value  of  about  0.111,  this  may  furnish  the 
necessary  indication. 

The  use  of  artificial  intelligence  techniques  to 
evaluate  the  residuals  and  trigger  all  of  the  manipulations 
of  the  moving  bank  has  been  suggested.  This  has  the 


advantage  of  looking  directly  at  the  real  source  of 


true  system  information.  However,  artificial  intelligence 
techniques  tend  to  be  computationally  intensive,  which  is 
the  opposite  of  the  goal  of  moving-bank  MMAE  as  opposed  to 
fnll-bank  NNAE,  so  there  is  a  tradeoff  inherent  in  this 
proposal . 

Finally,  issues  which  were  identified  prior  to  this 
thesis  still  need  to  be  addressed.  These  include  the 
evaluation  of  the  controllers  for  a  realistic  system 
(which  had  to  be  abandoned  here).  Also  the  study  of  slowly 
varying  uncertain  parameter  values,  and  uncertain  parameter 
values  not  equal  to  the  discrete  parameter  values  was  not 
addressed  here. 

VI. 4.  Summary 

This  thesis  has  addressed  the  evalnation  of  moving-bank 
multiple  model  adaptive  estimation  and  control  algorithms, 
as  applied  to  a  model  of  a  large  flexible  spacecraft.  It 
was  seen  that,  although  the  model  turned  out  not  to  require 
adaptive  control,  that  moving-bank  MHAE  could  still  provide 
estimates  of  the  uncertain  parameters.  Based  on  the 
performance  of  the  moving  bank  estimator,  several 
refinements  in  the  move  and  contraction  decision  logics  were 
made.  In  addition,  it  was  revealed  that  ambiguity  function 
analysis  can  be  an  invaluable  tool  in  designing  a  moving 
bank  estimator.  The  moving  bank  technique  continues  to  show 
great  promise  as  a  practical  adaptive  e s t i m a t i on/ c on t r o 1 
design  strategy,  and  research  in  this  area  should  continue. 


Appendix  A 


Detailed  Model  Deve lopment 

The  equations  for  the  elements  of  the  mass  and 
stiffness  matrices  were  given  in  Chapter  III  in  integral 
form.  They  are  repeated  here  for  clarity,  then  the 
simplified  closed  form  is  presented.  The  equations  for  the 
control  and  Measurement  matrices  are  also  developed  in  a 
more  complete  form.  Then  samples  of  the  M  (mass),  K 
(stiffness),  B  (control),  and  H  (measurement)  matrices  are 
given,  followed  by  two  tables  showing  how  the  free  system 
eigenvalues  of  the  second  and  fourth  bending  modes  vary 
throughout  the  parameter  space. 

The  equations  for  the  elements  of  the  mass  matrix  are 
repeated  from  Chpter  III  for  i  =  l,...,p  ;  and  j  =  l,...,p  : 

M(l.l)  =  2  (Ic  +  I2c) 

R+L 

M(l,i+1)  =  2pth  J  x0i  dx  +  2m2 ( R+L) 0 ^ ( z=x-R=L) 

R 

+  2I2c0'i(z=x-R=L) 

R+L 

M(I,i+p+l)  =  2pth  1  y0£  dy  +  2m2 ( R+L) 0 ^ ( x=y-R=L) 

R 

+  2I2c0'i(z=y-R=L) 

M ( i  +  1,1)  =  M ( 1 ,  i  +  1 ) 

M( i+p+1,1)  =  M( 1 , i+p+1 ) 
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M( i  +  1 , j  +1 )  =  2p  th  J  0.0.  dx 

J 

R 

+  2«2#i(z=i-R=L)»j (z=x-R=L) 

R+L 

M( i+p+1 , j +p+l )  =  2pth  J  0^0j  dy 

R 

+  2ai2^  i  (  z=y-R=L)  0  j  (  z=y-R=L) 

+  2I2c#'i(x=y-R=L)a'j  (z  =  y-R=L) 

M( i+1 , j +p+l )  =  0 
M(i+l+p,i+l)  =  0 

where  0^  is  a  function  describing  the  bending  modes  of  the 
arms ,  de f ined  a s : 

0^(2)  =  1  -  cos(  inz/L)  +  ( 1 /2 ) ( -1 )  i  +  2 ( in z /L) 2 
z=x-R(ory-R) 

and  J' j  is  the  first  derivative  of  0^  with  respect  to  z. 

The  pysical  dimensions  of  the  model  t,  h,  mj .  I 2C»  Ico» 
R,  and  L,  are  as  shown  in  Figure  III-l.  In  addition  p  is 
the  mass  density  of  the  arms,  and: 

Ic  =  ( 1 / 2 ) I Q  +  2  I j  +  2m2(R+L)2 

R+L 

Ij  =  pth  J  x2  dx 

R 
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The  elements  of  the  stiffness  matrix  are  given  for 
i  =  ;  and  j  =  ;  by: 

1(1.1)  -  0 
K( i+1.1)  =  0 
K ( i+l+p , 1 )  =  0 
K( 1 , j+1)  =  0 
K(l,j+p+l)  *  0 

L 

K(i+l,j+l)  =  2EI  J  #i"(z>0j"(s)  dz 

0 

K ( i+l+p , j +l+p )  =  K(i+l,j+l) 

K( i+l+p, j+1)  *  0 

K(i+l,j+l+p)  =  0 

where  E  is  the  Modulus  of  Elasticity  of  the  arms,  and  I  is 
the  area  moment  of  inertia  based  only  on  t  and  h  as: 

I  =  (l/12)ht3 

and  is  the  second  derivative  of  with  respect  to  z. 

When  the  necessary  substitutions  are  made  and  the 
integrals  are  evaluated,  the  simplified  equations  become: 

M ( 1 , 1 )  =  I„  +  4th(L3/3  +  RL2  +  R2L) 

+  m2(R2  +  2RL  +  L2)  +  2  I2(J 

M ( 1 , i  +  1)  =  2  p  t  h  {L2 / 2  +  [(-l)(i  +  1)  +  1 ] L/ ( in) 2 

+  (-l)(i+1)(in)2L(L/8  +  R/6)  +  RL) 

+  2m2(R  +  L) [ 1  +  (-l)(i+1)(l  +  (in)2/2)] 

+  2I2c(-1) ( i+1) I2n4/L 

M( 1 ,  i  +  p  +  1 )  =  M(  i  +  p+1 , 1 )  =  M( i+1  , 1 )  =  M ( 1 ,  i  + 1 ) 


M(i+l,j+l)  =  2p  th  (L+[(-l) ( i  +  1)  i2  +  (-1) ( J+1) j2n2L/6 

+  (-l)(i+J)L[  Ri2/Ri2  +  Ri2/Rj2 
+  (ijn)2/20]J 

+  2b2[1  +  (-l)(i  +  1)(l  +  (in)2/2)] 

[1  +  (-l)(i  +  D(i  +  (  j  n)  2  /  2  )  ] 

+  2I2c(-l)(i+j)(ijn/L)2 
+  {0:  i/j;  pthL:  i=j} 

M( i+p+1 , j +p+l )  =  M(i+l,j+l) 

M(i+l,j+p+l)  =  M(i+p+l,j+l)  =  0 

and 

K ( i  +  1 , j  +1 )  =  Et3h(-l)(i  +  j)L(ijn)2/(6L4) 

+  {0:  i#j f  Et3  h( ijn)2/ (12L3) :  i= j } 

K(l.l)  =  K  (  i  +  1 , 1 )  =  K  ( 1 ,  i  +  1 )  =  I(i+p+l,l)  «  K(l. i  +  p  +  1) 
=  K(i+l,j+p+l)  =  K(i+p+l,j+l)  =  0 


The  B  matrix  is  given  by: 


B 


12  2 
0  2£'i(z=L/2)  0 

0  0  2£'i(z=L/2) 


where  £^(z)  is  of  dimension  p  (the  number  of  modes  of 
interest),  and 

0  »  A ( z  =  L/2 )  = 


(  in/L)  sin(  iit/2)  +  ( -1 )  (  i  +  1 )  (  in  )  2  /  8 


The  H  matrix  is  given  by: 


H  = 


1 

0 

0 

0 

0 


0T  0T 
4Ti(z=L/2)  0T 
ftTi(z=L)  0T 

0T  4Ti(z=L/2) 

0T  ftTi(z=L) 


1  column  p  columns  p  columns 
H  is  thns  a  5  x  N  matrix,  and 

0a(z  =  L/2)  =  1  -  co s ( in/2)  +  0 . 5 ( - 1 ) ( i  +  1 > (  in ) 2 

and 

0i(z  =  L)  =  1  -  (-1) ( i  +  1) (  in)2 


When  these  expressions  are  evaluated  for  the  constant 
values  given  in  Table  III-l  (this  equates  to  the  (6,6)  point 
in  the  p/E  parameter  space),  and  for  p  =  2,  the  following 
matrices  result: 
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0.0 
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0.0 

-9.87 

0.0 

0.0 

0.0 

4.04 
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2.23 
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• 

o 

4.93 

-17.73 

As  seen  in  Chapter  III  the  free  system  eigenvalues  of  the 
system  can  be  found  via  Equation  (III-2): 

Kx  =  X  j 

The  eigenvalues  of  the  second  and  forth  bending  modes 
(the  unison  bending  modes  to  be  controled  in  this 
application)  are  tabulated  in  tables  A- 1  and  A- 2  for  all  the 
discrete  values  of  p  and  E  in  the  parameter  space. 
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F igure  B-14  . 


(  cont .  ) 


Appendix  C. 
SifislAlifiA  SflXi.wj.je. 


Four  main  problems  were  used  for  this  inve s i ga t ion.  A 
preprocessor  (PROGRAM  SETUPS),  used  to  create  the  parameter 
space;  a  primary  processor  (PROGRAM  BANE),  used  to 
generate  the  simulation  data;  a  post  processor  (PROGRAM 
RESLT),  used  to  plot  the  results;  and  a  program  to  perform 
ambiguity  function  analysis  (PROGRAM  AMBIG).  For  each 
program  this  Appendix  contains  a  structure  chart,  and  the 
program  and  subroutine  headers  (which  describe  each  module 
and  present  pseudocode  for  the  algorithms  performed).  The 
actual  FORTRAN  code  is  not  included  but  is  available  through 
Dr.  P.  S.  Maybeck  at  the  Air  Force  Institute  of  Technology, 
Department  of  Electrical  Engineering.  The  programs  make  use 
of  four  Libraries  available  on  the  Aeronautical  Systems 
Division  CDC  Cyber  computer:  IMSL5,  LQGLIB,  DISSPLA,  and 
CCLIB. 
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DESCRIPTION:  THIS  PROGRAM  CALCULATES  ALL  OP  TIE  CONSTANT  MATRICES 
REQUIRED  BT  PR06RAM  BANK  AND  AMBI6. 

AUTBOR  :  KARL  HENTZ  AND  PAUL  0.  FILIOS 
VERSION:  3.03 
DATE  :  13  AU6  S3 

INPUT:  FILE  aPARANTRa  -  CONTAINS  TIE  DISCRETE  PARAMETER  VALUES 

FOR  TIE  TWO  UNCERTAIN  PARAMETERS 

OUTPUT:  FILE  aSPACEa  -  CONTAINS  ALL  OF  TIE  CONSTANT  MATRICES  FOR 

EACB  DISCRETE  PARAMETER  POINT 


T  I  I  S 


TBE  MAIN  PRR06RAM 


PSEUDOCODE: 

READ  IN  TIE  PARAMETER  VALUES 
FOR  EACB  PARAMETER  POINT 

COMPUTE  TIE  CONTINUOUS  TIME  SYSTEM  MATRICES 
DISCRETIZE  TIE  NOISE  INPUT  MATRIX 
DISCRETIZE  TIE  PLANT  AND  INPUT  MATRICES 

COMPUTE  TBE  COVARIANCE  MATRIX  AND  TBE  FILTER  6AIN  MATRIX 
COMPUTE  TBE  A  INVERSE  MATRIX 
COMPUTE  TIB  DETERMINANT  OF  A 
SAVE  TBE  MATRICES  TO  LAI6BR  STORAGE  ARRAYS 
CONTINUE 

WRITE  THE  MATRICES  TO  TIE  OUTPUT  FILE 
END 


SEBBOETINE  ST STEM  (DEN.  E.  FMAT,  BEAT,  BEAT) 


OBBOUTINE  SYSTEM 


DESCBIPTION : 

THIS  PEOOBAM  COMPETES  M  (MASS).  E  (STIFFNESS).  SVBB  (CONTIOL) 
AND  SOBH  (MEASDBBMENT)  MATBICES  FOB  TBE  DEAPEE/EPL  SIMPLIFIED 
LAM6E  SPACE  STBECTEBE.  TBEN  OSES  TBEM  TO  COMPETE  TBE  STANDAED 
STATE  SPACE  F.  B.  AND  B  MATBICES . 

AETBOB:  PAEL  0.  FILIOS 

VEBSION :  1.0 

DATE:  29  JEL  1SSS 

INPET:  DEN.  E  -  TME  VALDES  FOB  TBE  TBO  DNCBBTAIN  PABAMETBBS 

ODTPET:  FMAT.  BMAT.  MMAT  -  CONTINOE8  TIME  MATBICBS  TMAT  DESCBIBE 

TBB  SYSTEM 


PSEUDOCODE: 

BEGIN 

CALCDLATE  TME  DIMENSIONS  OF  TBE  MATBICES 

CALL  SEBBOETINE  TMAT  COMPETES  TBE  M  AND  E  MATBICES 

CALL  EIGENVALEE/EIGENVECTOB  BODTINE 

FOBM  SMAT  MATE IX 

FOBM  SDBB  MATBIX 

FOBM  SDBM  MATBIX 

CALCDLATE  SMAT  TBAN8P0SE  •  SEBB  -  NBVB 
FOBM  BMAT  MATBIX 
CALCDLATE  SDBB  •  8  -  NBVB 
FOBM  BMAT  MATBIX 
END 


o  u  u 
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DESCRIPTION:  THIS  PIOGEAM  SIMULATES  A  MOVING  DANK  MULTIPLE  MODEL 
ADAPTIVE  KALMAN  PILTBE.  BASED  ON  TIE  TEUTI  MODEL  DB8CEIBBD  BT 
SUBEOUTINE  TEUE  AND  CONTBOLS  TIE  SYSTEM  BAUD  ON  TIE 
CONTEOL  SYSTEM  DESCRIBED  BY  SUBEOUTINE  CNTEL. 

AUTBOB  :  KABL  IBNTZ  AND  PAUL  6.  PII IOS 

VEBSION :  2.4 

DATE  :  02  OCT  SS 

INPUT:  PILE  ' SPACE  * p  VEITTEN  BY  TIE  SETUPS  PE06EAM.  CONTAINS  TIE 

MATEICES  POE  BACI  PILTEI  IN  TIE  PAEAMETBE 
SPACE. 

PILE  * PAIAMTE ' p  CONTAINS  TIE  DI8CEETE  PAIAMETEB  VALUES  POE 
EACB  VAEYING  PAEAMETEE. 

PILE  *CONDIT>  CONTAINS  TIE  INPUT  CONDITIONS  FOB  TIE 
SITUATION  UNDBE  STUDY. 

OUTPUT:  FILE  'OUTl’n  CONTAINS  TIE  IA*  DATA  FOE  TIE  EESLTS 

PI06BAM  TO  PUT  INTO  GEAPI  POEM. 

FILE  *OUT2 ' p  CONTAINS  TIE  DATA  GBNEEATED  BY  TIE  FIEST 
MONTE  CAELO  BUN  IN  A  POEM  SUITABLE  POE 
FEINTING. 

PILE  ' STATES >  CONTAINS  TIE  VALUES  OP  TIE  SYSTEM  STATE 

VAEIABLBS  AT  EACB  SAMPLB  TIME  OP  TBE  FUST 
MONTE  CAELO  BUN. 

NOTB:  THIS  PIOGEAM  BAS  BEEN  MODIFIED  TO  USE  ONLY  AN  APPBOXINATION 
POE  TIE  BANK  VBIGITIN6  PEOBABILITIES  DUE  TO  NUMEIICAL  DIFFICULTIES 
IN  COMPUTING  ADET  IN  THE  SETUPS  PIOGEAM.  ONLY  SUBEOUTINE  aFii 
IS  AFFECTED. 


PSEUDOCODE: 

READ  THE  INPUT  VALUES 
INITIALIZE  OUTPUT  FILES 
FOR  EACR  MONTE  CARLO  RUN 

INITIALIZE  TBE  TRUE  PARAMETERS 
INITIALIZE  TBE  BANE  VARIABLES  AND  MATRICES 
FOR  EACB  SAMPLE  PERIOD 

PROPAGATE  TBE  TRUE  SYSTEM 
PROPAGATE  TBE  ESTIMATE 
TAKE  A  MEASUREMENT 
UPDATE  TBB  BANK  ESTIMATES 
DETERMINE  TBE  BANK  PROBABILITIES 
DETERMINE  TEE  NEB  E8TIMATE 
COMPUTE  TBE  CONTROL  INPUTS 
FOR  TBE  FIRST  RUN 

OUTPUT  TBE  PROBS,  STATES.  AND  TRRESBOLD  VARIABLES 
MOVE.  CONTRACT.  OR  EXPAND  TBE  BANK  IF  NBCESSART 
COMPUTE  TBE  PERFORMANCE  STATISTICS 
OUTPUT  TBE  STATISTICS  OF  INTEREST 
LOOP 
LOOP 
END 


SUBROUTINE  CNTRL  ( AHI ST , XBATA , U , TIM . ALPS AB , XBATB , PIOB ) 


SUBROUTINE  CNTXL 


DESCRIPTION :  THIS  ROUTINE  CONFUTES  TBE  CONTROL  INPUT 

AUTHOR  :  CARL  HENTZ  AND  PAUL  G.  FILIOS 

VERSION  :  3.0 

DATE  :  17  SEP  S3 

PARAMETERS  PASSED: 

AHI ST:  THE  ESTIMATES  OP  THE  PARAMETER  VALUES  POR  THE  LAST  10 
TIME  PERIODS.  USED  POR  ADAPTIVE  CONTROL 
XHATA:  THE  WEIGHTED  AVERAGE  OP  STATE  ESTIMATES.  USED  FOR  ADAPTIVE 
CONTROL 

U:  THE  CONTROL  INPUT  VECTOR.  COMPUTED  HERE 

TIM:  THE  CURRENT  SIMULATION  TIMB  PERIOD 

ALPHAS:  THE  INDIVIDUAL  PARAMETER  ESTIMATES  OP  THE  BANC  PILTERS 
XHATB :  THE  INDIVIDUAL  STATS  ESTIMATES  OP  THE  BANC  FILTERS 

PROB:  THE  PROBABILITT  WEIGHTINGS  OP  THE  BANC  FILTERS 


PSEUDOCODE: 

•  CASE  OP 

•  DITHER  SIGNAL 

•  FIXED  GAIN  CONTROLLER 

•  COMPUTE  CONTROL 

•  VARIABLE  GAIN  CONTROLLER 

•  SELECT  GAIN 

•  COMPUTE  CONTROL 

•  MMAC 

•  FOR  EACH  FILTER 

•  COMPUTE  CONTROL 

•  ACCUMULATE  WEIGHTED  AVERAGE 

•  LOOP 

•  END 
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SUBROUTINE  TRUE  (X . U . ATRUE, RNSTP.TIM) 


SUBROUTINE  TRUE 


DESCRIPTION:  THIS  ROUTINE  PROPAGATES  THE  TRUE  SYSTEM 
FROM  ONE  SAMPLE  TIME  TO  TRE  NEXT.  ADJUSTING  TEE  TRUE 
PARAMETER  POINT  VBEN  DESIRED 

AUTHOR  :  EARL  RENTZ  AND  PAUL  G.  FILIOS 

VERSION  :  4.0 

DATE  :  23  SEP  S3 

PARAMETERS  PASSED: 

X:  THE  TRUE  SYSTEM  STATES 

U:  TRE  CONTROL  INPUT  VECTOR 

ATRUE:  TRE  TRUE  PARAMETER  POINT 

RN8TP :  TIME  OP  TRE  END  OF  TRE  SIMULATION 

TIM:  TRE  CURRENT  SIMMULATION  TIME  PERIOD 


PSEUDOCODE: 

IF  TME  TRUE  PARAMETER  VARIES 

GET  TRE  NEV  MATRICES  FROM  TRE  PARAMETER  SPACE 
MULTIPLY  THE  STATES  BY  THE  STATE  TRANSITION  MATRIX 
MULTIPLY  THE  CONTROL  BY  TRE  TRUE  RD  MATRIX 
GET  A  NOISE  VECTOR 
MULTIPLY  TRE  NOISE  VECTOR  BY  RD 
ADD  ALL  THE  PRODUCTS  FOR  THE  NEV  TRUE  STATES 
END 
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SUBROUTINE  IV (V) 


SUBROUTINE  RV 


DESCRIPTION:  GENERATES  A  WHITE  6AUSIAN  RANDOM  VECTOR. 

AUTHOR  :  EARL  HENTZ  AND  PAUL  G.  FILIOS 

VERSION:  2.0 

DATE  :  2 <  JUN  85 

PARAMETERS  PASSED: 

V:  A  RANDOM  NOISE  VECTOR 


PSEUDOCODE; 

FOR  EACH  ELBMENT  OF  THE  NOISE  VECTOR 
DO  12  TIMES 

ACCUMULATE  A  UNIFORMLY  DISTRIBUTED  RANDOM  NUMBBR 
LOOP 

SUBTRACT  SI1 
LOOP 
END 


SUBROUTINE  PBOPOT ( XIATB ,  PUB,  BDB .  U > 


SUBROUTINE  PBOP6T 


DESCBIPTION :  TBI3  ROUTINE  PROPAGATES  TBE  STATE  ESTIMATES  OF  EACE 
FILTER  IN  TBE  MOVING  BAN!  FROM  TBE  END  OF  ONE  SAMPLE  PERIOD  TO 
JUST  BEFORE  TBE  MEASUREMENTS  ARE  TAKEN  OF  TIE  NEXT. 

AUTBOI  :  KARL  BENTZ  AND  PAUL  G.  FILIOS 
VERSION:  2.2 
DATE  :  3  AUG  S3 

PARAMETERS  PASSED: 

XBATB :  CURRBNT  STATE  ESTIMATE  FOR  EACB  FILTBR  IN  TIE  BANK.  TBESE 
VALUES  ABE  ALTERED 

PUB:  PBI  MATRIX  FOR  EACB  FILTER  IN  TIB  BANK 

BDB:  DISCRETE  INPUT  MATRIX  FOR  EACB  FILTER  IN  TIE  BANK 

U:  CONTROL  INPUT  VECTOR 


PSEUDOCODE: 

FOR  EACB  FILTER 

MULTIPLY  PBI  BT  XIAT 
MULTIPLY  BD  BY  U 
ADD  TIE  PRODUCTS  TOGETIER 
LOOP 
END 
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SUBROUTINE  MBAS<X. HMATT, Z.  R) 


SUBROUTINE  NBAS 


DESCRIPTION:  TBIS  ROUTINE  TAKES  A  MEASUREMENT  OF  TME  TRUE  STSTEM. 


AUTMOE  :  KARL  MENTZ  AND  PAUL  0.  F1LIOS 
VERSION:  3.02 
DATE  :  3  AUS  IS 


PARAMETERS  PASSED: 

X :  TRUE  STATES 

MNATT :  TRUE  MEASUREMENT  MATRIX 

Z:  TME  MEASUREMENT  RETURNED 

R:  TME  COVARIANCE  OF  TEE  MEASUREMENT  NOISE 


PSEUDOCODE: 

CALL  THE  P8UED0- RANDOM  NUMBER  0ENERATOR  FOR  THE  V  VECTOR 
MULTIPLE  TEE  V  VECTOR  BT  THE  SQUARE  ROOT  OF  R 
MULTIPLY  X  BT  MMATT 
ADD  TME  NOISE  TO  THE  MEASUREMENT 
END 


<S£vt  v  VivOvYv.v  ivi*£>£>X  ■$  >$■ 


SUBROUTINE  EES ID ( Z , ZIATB , MNATB . EESB , LIKE . AINVB ) 


SUBROUTINE  RESID 


DESCRIPTION:  THIS  ROUTINE  COMPUTES  TBE  RESIDUALS  FOR  BACH  OF  THE 
ELEMENTAL  FILTERS  IN  TBE  SLIDING  BANE.  AND  TEE  LIKELIHOOD 
QUOTIENT  FOR  EACH  FILTER 

AUTBOR  :  KARL  HBNTZ  AND  PAUL  6.  FILIOS 
VERSION:  2.03 

DATE  :  3  AUO  S3 

PARAMETERS  PASSED: 

Z:  TBE  MEASUREMENTS  TAKEN 

1HATB :  THE  CURRENT  STATE  ESTIMATES  FOR  EACH  FILTER  IN  THE  BANK 
HMATB :  THE  MEASUREMENT  MATRICES  FOR  BACB  FILTER  IN  TBE  BANK 
RESB :  THE  RESIDUALS  RETURNED 

LIKE:  THE  MINIMUM  LIKELIHOOD  QUOTIENT  RBTURNED 
AINVB:  THB  A  INVERSE  MATRIX  FOR  EACH  FILTBR  IN  THE  BANK 


PSEUDOCODE: 

FOR  EACH  FILTER 

MULTIPLY  HMAT  BT  XHAT  TO  DETERMINE  Z  EXPECTED 
COMPUTE  THE  RESIDUAL  -  Z  -  Z  EXPECTED 
FORM  THE  QUADRATIC  R  •  AINV  •  R(TRANSPOSE) 

IF  THAT  IS  LESS  THAN  THE  MINIMUM  LIKELIHOOD  QUOTIENT  THEN 
MAKE  IT  THE  NB*  MINIUM  LIKELIHOOD  QUOTIENT 

LOOP 

END 
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SUBROUTINE  UPDTE ( RESB , XHATB .PMINB , AINVB , CENATB . TIM ) 


SUBROUTINE  UPDTE 


DESCRIPTION:  TBIS  ROUTINE  UPDATES  EACB  OF  TIE  FILTERS  IN  TIE 
SLIDING  BANK  BITS  TIE  RESULTS  OF  TIE  MEASUIEIINTS. 

AUTBOI  :  KARL  BENTZ  AND  PAUL  G.  FILIOS 
VERSION:  2.03 
DATE  :  3  AUG  S3 

PARAMETARS  PASSED: 

RESB:  TIE  RESIDUALS  FOR  EACB  FILTER  IN  TIE  BANK 

XIATB :  TIE  CURRENT  STATE  ESTIMATES  FOB  EACB  FILTER  IN  TIE  BANK 

PMINB:  TIE  COVABIANCE  MATRIX  BBFOIB  UPDATE  FOR  EACB  FILTER 

AINVB:  TBE  A  INVERSE  MATRIX  FOR  EACI  FILTER  IN  TIE  BANK 

CKMATB :  TIE  KALMAN  FILTER  GAIN  MATRIX  FOR  BACB  FILTER  IN  TIB  BANK 

TIM:  TBE  CURRENT  TIME  IN  TBE  SIMULATION 


PSEUDOCODE: 

FOR  EACB  FILTER 

MULTIPLE  TBE  GAIN  MATRIX  BT  TBE  RESIDUALS  TO  GET  TIE  C RANGE  IN 
TBE  STATE  ESTIMATES 

ADD  THE  CHANGE  TO  TBE  CURRENT  STATE  ESTIMATES 
LOOP 
END 


14  7 


■  s'  VuV 


SUBROUTINE  PROBL ( RESB . AI NVB , ADETB . PHOB , PLOW ) 


•••»*•••••••••••••••••••••••••• 

SUBROUTINE  PROBL 


DESCRIPTION:  TRIS  ROUTINE  COMPUTES  TEE  PROBABILITY  VEIORTING 
FACTORS  FOR  EACH  FILTER  IN  TIE  BANE 

AUTHOR  :  EARL  BENTZ  AND  PAUL  6.  FILIOS 
VERSION:  2.00 
DATE  :  20  JUL  IS 

PARAMETERS  PASSED: 

RESB:  THE  RESIDUALS  FOR  EACH  FILTER  IN  THE  BANE 

AINVB :  THE  A  INVERSE  MATRIX  FOR  EACH  FILTER  IN  THE  BANE 

ADETB:  THE  DETERMINATE  OF  TEE  A  MATRIX  FOR  BACH  FILTER 

PROB:  THB  CUREENT  PBOBABILITY  VBIGHTINOS  FOR  EACH  FILTER.  RETURNED 

PLOR:  THE  PROBABILITY  WEIGHTING  LOVER  BOUND 


PSEUDOCCDE: 

FOR  EACH  FILTER 

GET  THE  PROBABILITY  DENSITY  FUNCTION  FOB  EACH  FILTER 
MULTIPLY  IT  BY  THE  LAST  PROBABILITY  WEIGHTING 
ACCUMULATE  THE  TOTAL  VALUES  FOB  USE  IN  THE  DNOMINATOR 
LOOP 

DIVIDE  ALL  PROBABILITY  VALUES  BY  THE  ACCUMULATED  TOTAL 
END 


SUBROUTI NE  F ( RES . AIHV. ADET , DNSTT ) 


SUBROUTINE  F 


DESCRIPTION:  TBIS  ROUTINE  COMPUTES  TEE  PROBABILITY  DENSITY 
USED  IN  C0MPUTIN6  TEE  BANE  FILTER  PROBABILITIES. 

AUTBOR  :  KARL  BENTZ  AND  PAUL  G.  FILIOS 
VERSION:  2.00 
DATE  :  20  JUL  SS 

PARAMETERS  PASSED: 

RES  :  INPUT  RESIDUALS  FOR  ONE  FILTER 

AINV:  INPUT.  TRE  INVBRSE  A  MATRIX  FOR  ONE  FILTER 

ADET:  INPUT.  TBB  DETERMINANT  OF  TRE  A  MATRIX  FOX  ONE  FILTER 

DNSTT:  OUTPUT.  TBB  PROBABILITY  DENSITY  COMPUTED 

NOTE:  BECAUSE  OF  NUMERICAL  DIFFICULTIES  IN  CALCULATING  ADET 
DNSTY  IS  INCOMPLETELY  CALCULATED.  TEE  TRUE  VALUE  OF  DNSTT 
CAN  BE  IMPLEMENTED  BY  REMOVING  TRE  *•«  COMMENT  SYMBOL  FROM 
TRE  LINE  BEFORE  THE  RBTURN  STATEMENT. 


PSEUDOCODE: 

MULTIPLY  THE  RESIDUALS  TRANSPOSE  BY  A  INVERSE 
MULTIPLY  THE  PRODUCT  BY  TBE  RESIDUALS 
DIVIDE  THE  PRODUCT  BY  -2 
RAISE  E  TO  TBE  RESULTING  QUOTIBNT 

SCALE  BY  2*PI  ••  S  •  SQUARE  ROOT  OF  TBE  DETERMINANT  OF  A 
END 


SUBROUTINE  A VO ( XBATB . ALPBAB .PMINB , PALPB, PBOB . XBATA . ABI ST) 


»  SUBBOUTINE  AVC 


DESCRIPTION :  TB1S  ROUTINE  TAXES  TBE  VEIGRTBD  AVEBACE  OF  TBE 
INDIVIDUAL  FILTEB’S  STATE  AND  PABANBTEB  ESTINATES 

AUTBOR  :  KARL  BENTZ  AND  PAUL  0.  FILIOS 
VERSION:  2.00 
DATE  :  20  JOL  S3 

PARAMETERS  PASSED: 

XBATB:  TBE  STATE  ESTIMATES  FROM  TBE  FILTERS  IN  TBE  BANK 

ALPBAB:  TBE  PARAMETER  ESTIMATES  FROM  TBE  FILTERS  IN  TBE  BANK 

PMINB:  TBE  COVARIANCE  OF  TBE  PARAMETER  ESTIMATES  BEFORE  UPDATE 

PALPB:  TBE  COMPUTED  COVARIANCE  OF  TBE  OVERALL  PARAMETER  ESTIMATE 

RPROB :  TBE  PROBABILITY  VEIOBTINOS  FOR  BACB  FILTER 

XBATA:  TBE  AVERAOED  STATE  ESTIMATE  COMPUTED  AND  RETURNED 

ABIST:  TBE  BISTORT  OF  TBE  LAST  10  PARAMETER  ESTIMATES,  UPDATED 

AND  RETURNED 


PSEUDOCODE: 

SBIFT  TBE  OLD  VALUBS  OF  ALPBA  IN  TBE  BISTORT  ARRAT 
FOR  EACR  FILTER 

ACCUMULATE  XBAT  •  TBE  PBOBABILITT  VEIGBTING 
ACCUMULATE  ALPBA  •  TBE  PBOBABILITT  VEIGBTING 
LOOP 

COMPOTE  PALPB 
END 
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SOBBODTINE  DSEID(ABIST, PALPI, LIKE, A I NVB , PMINB , PBIB , BOB , XBATB , 
BMATB , ADETB , ALPBAB , CKBATB . PBOB . Z , ZBATA . I CENT , I  SIZE) 


SOBBODTINE  DSEIO 


DESCBIPTION :  TBIS  BODTINE  EXAMINES  TBE  SLIDING  BANE  AND  DESIDES 
IF  IT  NEEDS  TO  BE  MOVED,  EXPANDED,  OB  CONTBACTBD,  TBEN  CALLS 
TBE  APPBOPBIATE  BODTINE  TO  ACCOMPLI SB  ANT  NBCBSSABT  ACTION. 

ADTBOB  :  PADL  0 .  FILIOS 
VEBSION :  3.00 
DATE  :  1«  ADO  S3 

PABAMETEBS  PASSED: 

AINVB,  PMINB.  PBIB,  BOB.  XBATB ,  BMATB ,  ADETB,  ALPBAB.  CEMATB, 
PBOB:  TBE  ABBATS  TBAT  DBSCBIBE  TBE  FILTBBS  IN  TBE  BANE,  IF  TBE 
BANE  MOVES,  CONTBACTS.  OB  EXPANDS  TBESE  ABBATS  BILL  BE  DPDATED 
TO  BEFLECT  TBE  NEB  FILTEBS  IN  TBE  BANE 
PALPB,  LIEE,  PBOB:  TBE  VALDES  OF  TBESE  PABAMETEBS  ABE  DSED  IN  TBE 
DECISION  LOGIC 

I CENT.  I SIZE:  DESCBIBE  TBE  COBBENT  CBNTEB  AND  SIZE  OF  TEE  BANE 
AEIST;  CONTAINS  TBE  LAST  ESTIMATE  OF  TBE  DNCEBTAIN  PABAMETEBS, 
DSED  TO  DETBBMINE  TEE  NEB  BANE  CENTEB  AFTEB  A  CONTBACTION 
XBATA :  TBE  CDBBENT  STATE  ESTIMATE,  DSED  TO  INITIALIZE  NEB  FILTEBS 


PSEDDOCODE: 

IF  TBE  BANE  NEEDS  TO  CONTBACT 
CONTBACT 

IF  TBE  BANE  NEEDS  TO  MOVE 
MOVE 

IF  TBE  BANE  NEEDS  TO  EXPAND 
EXPAND 

END 


131 


SUBROUTINE  CNTRCT  < I  CENT , I  SIZE . PHIB , BDB , PMINB . AINVB , ADETB , CKMATB , 
PROB . BHATB . IBATB , XIATA . ALPHAS ) 


SUBROUTINE  CNTRCT 


DESCRIPTION:  THIS  ROUTINE  CONTRACTS  THE  SLIDING  BANK  TO  A  FINER 
DISCRETIZATION. 

AUTHOR  :  PAUL  G.  FILIOS 
VERSION:  1.00 
DATE  :  20  AUG  SS 

PARAMETERS  PASSED: 

ICENT.  ISIZE:  THE  NE*  CENTER  AND  SIZE  OF  THE  BANE 
PHIB . BDB . PMINB . AINVB . ADETB . CKMATB . PROB , HMATB . XHATB . ALPHAS :  THE 
ARRATS  DESCRIBING  THE  FILTERS  IN  THE  BANK 
XHATA.  TBE  CURRENT  STATE  ESTIMATE.  USED  TO  INITIALIZE  THE  NEW 
FILTERS  IN  THE  BANK 


PSEUDOCODE: 

ENSURE  ICENT  HILL  KBEP  SLIDING  BANK  INSIDE  PARAMETER  SPACE 
IF  NOT  ADJUST  ICENT 
REAS1GN  THE  PHI  MATRICES 
REAS1GN  THE  K  MATRICES 
REASIGN  THE  HD  MATRICES 
REASIGN  THE  H  MATRICES 
REASIGN  THE  P  MINUS  MATRICES 
REASIGN  THE  AINV  MATRICES 
REASIGN  THE  A  DETERMINATES 
REASIGN  THE  ALPHAS 
REASIGN  TBE  PROBABILITIES 
ASIGN  XBAT  TO  NER  FILTERS 
END 


SUBROUTINE  EXPAND ( I CENT, I SIZE , PUB . BOB , PHINB , AINVB , ADETB , CEMATB , 
*  PROB , HMATB . XIATB .XHATA, ALPHAS) 


SUBROUTINE  EXPAND 


DESCRIPTION:  THIS  ROUTINE  EXPANDS  THE  SLIDING  BANK  TO  ITS 
LARGEST  DISCRETIZATION. 

AUTBOR  :  PAUL  G.  FILIOS 
VERSION:  1.00 
DATE  :  20  AUG  S5 

PARAMETERS  PASSED: 

ICBNT.  I SIZE :  TBE  NE*  CENTER  AND  8IZE  OP  THE  BANK 
PHIB, BDB , PMINB .AINVB . ADETB . CKMATB . PROB . HMATB , XHATB , ALPHAS :  THE 
ARRATS  DESCRIBING  THE  FILTERS  IN  TBE  BANK 
XHATA.  THE  CURRENT  STATE  ESTIMATE.  USED  TO  INITIALIZE  THE  NER 
FILTERS  IN  THE  BANK 


PSEUDOCODE: 

SET  ISIZE  TO  4 

SET  THE  BANK  CENTER  TO  THE  PARAMETER  SPACE  CENTER 
RESET  PHIB.  BDB,  HMAT,  CKMAT,  PMINB.  AINVB.  ADETB. 
ALPHAB.  PROB.  AND  XHATB 

END 
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SUBROUTINE  MOVE  < ICENT, I  SIZE , PHIB , BOB . PMINB, AINVB , ADETB. CKMATB, 
♦  PROS , BMATB . ZBATB , ZBATA . ALPBAB , NOVDIR , Z ) 


SUBROUTINE  HOVE 


DESCRIPTION:  TBIS  ROUTINE  MOVES  TBE  SLIDDNG  BANX  ONE  8TBP  IN  ANT 
DIRECTION 

AUTHOR  :  PAUL  G.  FILIOS 
VERSION:  3.00 
DATE  :  20  AUG  S3 

PARAMETERS  PASSED: 

ICENT.  ISIZE :  THE  NEB  CENTER  AND  SIZE  OP  TBE  BANK 
PHIB . BDB .PMINB .AINVB , ADETB , CKMATB , PROB . BMATB . XBATB . ALPBAB :  TBE 
ARRAYS  DESCRIBING  THE  FILTERS  IN  TBE  BANK 
ZBATA,  TBE  CURRBNT  STATE  ESTIMATE.  USBD  TO  INITIALIZE  TBE  NEV 
FILTERS  IN  TBE  BANK 


PSEUDOCODE: 

COMPUTE  NEV  BANK  CENTER 

IF  MOVE  PUTS  PART  OF  BANK  OUTSIDE  PARAMETER  SPACE 
THEN  RETURN 

RESET  PHIB.  BDB.  BMAT,  CKMAT,  PMINB.  AINVB,  ADETB 
ALPHAS,  AND  ZBATB 
ASSIGN  NEV  PROBABILITIES 
END 
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PROGRAM  RESLT 


PROGRAM  RESLT 


DESCRIPTION:  THIS  PROGRAM  IS  TIE  POST  PROCESSOR  FOR  PROGRAM  BANK. 
ORGANIZING  THE  RA*  DATA  INTO  A  FORM  THAT  CAN  BE  PLOTTED  BT  THE 
CALCOMP  PLOTTER. 

AOTHOR  :  KARL  HENTZ  AND  PAUL  G.  FXLIOS 

VERSION:  2.3 

DATE  :  20  OCT  S3 

INPOT:  FILE  '0UT1'  CREATED  BT  PROGRAM  BANK. 

OUTPUT:  FILE  ’TAPB99  *  READ!  FOR  ROUTING  TO  THE  CALCOMP  PLOTTER. 


THIS  IS  THE  MAIN  PROGRAM 


•  PSEUDOCODE: 

•  INITIALIZE  INPUT  ARRAYS  TO  0 

•  READ  IN  EAR  DATA.  SUMMING  THE  MONTE  CARLO  RUNS 

•  NORMALIZE  DATA  BT  THE  NUMBER  OF  MONTE  CARLO  RUNS 

•  SET  UP  THE  PLOT  TITLES 

•  PUT  THE  DATA  INTO  THE  PLOT  ARRATS 

•  PLOT  THE  DATA 

•  END 
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SUBROUTINE  SLOTS  < DATA, NPNTS . LTITL. TTITL, MODE . NN. X. T. NS. IU . 
♦  MOVE , EXPO, CONTHB . NBR . IPVAB . PSTIT) 


SUBROUTINE  SLOTS 


DESCRIPTION:  TBIS  ROUTINE  PLOTS  THE  DATA  PREPARED  IN  THE  MAIN 
PROGRAM  USING  CALCOMP  PLOTING  ROUTINES. 

AUTHOR  :  EARL  HENTZ  AND  PAUL  6.  FILIOS 

VERSION:  2.2 

DATE  :  21  OCT  SS 

PARAMETERS  PASSED:  (ALL  ARE  INPUT) 

DATA:  TME  DATA  POINTS  TO  BE  PLOTED 

NPNTS:  TEE  NUMBER  OP  DATA  POINTS 
LTITL:  THE  TITLE  OP  THE  PLOT 

TTITL:  AN  ARRAT  OF  TITLES.  ONE  FOR  THE  T  AXIS  OF  BACM  SUB-PLOT 

MODE:  THE  MODE  THE  8LID1N6  BANK  HAS  OPERATING  IN 

MOVE. 

EXPO, 

CONTHR:  THE  THRESHOLDS  USED  FOR  MOVE.  EXPAND,  AND  CONTRACT  DECISIONS. 
NBR:  THE  NUMBER  OF  MONTE  CARLO  RUNS  MADE. 

IPVAR:  CODB  SIGNIFYING  THE  VARIABLB  OF  INTEREST  BEING  PLOTTED. 

PSTART:  THE  TIME  IN  THE  MONTE  CARLO  RUNS  WHEN  THB  PLOT  STARTS. 


PSEUDOCODE: 

WRITE  TITLE  AND  NOTES 
FOR  EACH  SUB-PLOT 

DRAW  AND  LABLE  THE  AXIS 

DRAW  THE  MEAN  LINE 

DRAW  THE  PLUS  ONE  SIGMA  LINE 

DRAW  THE  MINUS  ONE  SIGMA  LINE 

ADJUST  STARTING  POINT  FOR  NEXT  SUB-PLOT 

END 
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PBOGBAM  ANBI6 


P  B  0  6  B  A  M  A  N  B  I  G 


OESCIIPTION:  PEKFOBNS  ABIGUITT  FUNCTION  ANALYSIS  OF  PABABBTBB  POINTS 
IN  TBE  PABAMETEK  SPACE  OF  A  SLIDING  BANE  MULTIPLE  MODEL  ADAPTIVE 
ESTIMATOB. 

AUTBOB  :  EABL  BBNTZ  AND  PAUL  G.  FILIOS 

VEBSION:  2.0 

DATE  :  20  ION  SS 

INPUT:  FILE  xAMSPACE*  CONTAINING  TEE  MATEICBB  DESCEIBING  TBB  PAEAMBTBB 
SPACE 

OUTPUT:  «TAPB9»«  PLOT  FILE  FOE  TIE  CALCOMP  PLOTTED 


T  M  I  8  IS  TBE  MAIN  PEOGBAM 


PSEUDOCODE: 

BEAD  IN  TBE  PAAMETEB  SPACE 
SET  UP  TBE  TEUTB  MODEL 
FOB  BACM  FILTBB 
SET  UP  TBE  FILTEE 
DETEBMINE  TBE  AMBIGUITY  FUNCTION 
LOOP 

PEEPABE  DATA  FOB  PLOTTING 
PLOT  DATA 
END 
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SOBBOUTINE  SETTB(P1 . P2 , PEIT, SDT. EMATT) 


SOBBOUTINE  SETT* 


DESCBIPTION :  8ETS  DP  TBE  TBDTE  MODEL  POE  TEE  PAEAMETEE  POINT  DESIEED 

ADTBOE  :  KAEL  BENTZ  AND  PAOL  6.  FILIOS 

VEE8ION :  2.0 

DATE  :  25  JDN  85 

PABAMETEBS  PASSED: 

PI ,P2 :  INDICES  TO  TIE  TIDE  PAEAMETEE  POINT  IN  TBE  PAEAMETEE  SPACE 

PUT:  TBDE  PB1  MATEIX,  EETDENED 

SDT:  TEDE  NOISE  INPUT  MATEIX.  EETDENED 

■MATT:  TEDE  MEASDEEMENT  MATEIX,  EETDENED 


PSEUDOCODE: 

SET  PUT  PEON  PAEAMETEE  SPACE 
OET  SDT  PBOM  PAEAMETEE  SPACE 
GET  EMATT  PEON  PAEAMETEE  SPACE 
END 


SOBBOOTI NE  SBTFI  ( PI .  P2  ,  HI .  HAT  ,  CHAT ,  PH  IN .  PPLOS .  AI NV ,  ABET ) 


SOBBOOTINE  SETF1 


DESCEIPTION :  SETS  OP  TIE  FILTEB  MODEL  FOE  TIE  PABAMETEB  POINT 
OESIEED. 

AOTIOB  :  KAIL  BENTZ  AND  PAOL  0.  FXLIOS 

VEBSION :  2.0 

DATE  :  25  JON  SS 


PABAMETEES  PASSED: 

PI.P2:  INDICES  TO  TIE  PAEAMETEE  POINT  DBSIEED 
PII :  STATE  TIAN81TI0N  MATIIX  FOE  DESIEED  POINT 
HAT:  MEASOBHBNT  MATEIX  FOE  DESIEED  POINT 
CHAT:  FILTBE  SAIN  MATEIX  FOE  DESIEED  POINT 
PMIN :  COVAEIANCE  OF  STATE  ESTIMATES  IEFOEB  OPDATES 

PPLOS:  COVAEIANCB  OF  STATE  ESTIMATES  AFTEE  OPDATES 
AINV:  A  MATEIX  FOE  DESIEED  POINT 
ADET:  DBTEEMINANT  OF  A  MATEIX 


PSEODOCODE: 

CET  PII  FIOM  PAEAMETEE  SPACE 
6ET  PMIN  FEOM  PAEAMETEE  8PACE 
GET  PPLOS  FEOM  PAEAMETEE  SPACE 
GET  CEMAT  FEOM  PABAMETEB  SPACE 
GET  IMAT  FEOM  PABAMETEB  SPACE 
GET  AINV  FEOM  PAEAMETEE  SPACE 
GET  ADET  FEOM  PABAMETEB  SPACE 
END 
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SUBROUTINE  !V(V) 


SUBROUTINE  RV 


DESCRIPTION:  GENERATES  A  f RITE  6AUSIAN  RANDOM  VECTOR. 

AUTBOR  :  EARL  BENTZ  AND  PAUL  6.  FILIOS 

VERSION:  2.0 

DATE  :  26  JON  <5 

PARAMETERS  PASSED: 

V:  A  RANDOM  NOISE  VECTOR 


PSEUDOCODE: 

FOR  EACB  ELEMENT  OF  TEE  NOISE  VECTOR 
DO  12  TIMES 

ACCUMULATE  A  UNIFORMLY  DISTRIBUTED  RANDOM  NUMBER 
LOOP 

SUBTRACT  SIX 
LOOP 
END 
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DESCRIPTION :  THIS  ROUTINE  TAKES  A  MEASUREMENT  OP  THE  TRUE  STSTEM 


ADTHOR  :  KARL  HENTZ  AND  PAUL  6.  FILIOS 
VERSION:  2.02 
DATE  :  3  AUG  SS 

PARAMETERS  PASSED: 

X:  TREE  STSTEM  STATES 

HNATT:  TRUE  STSTEM  MEASUREMENT  MATRIC 

Z:  MEASUREMENTS  TAKEN 

R:  COVARIANCE  OF  TEE  MEASUREMENT  NOISB 


PSEUDOCODE: 

GET  A  RANDOM  NOISE  VECTOR 

MULTIPLE  EACH  ELEMENT  BT  THE  SQUAREROOT  OF  R 

MULTIPLE  RMATT  TIMES  X 

ADD  THE  NOISE  TO  THE  PRODUCT 


SUBROUTINE  UPDTE ( Z .  IB .  MAT .  CHAT ) 


SUBROUTINE  UPDTE 


DESCRIPTION:  TBIS  ROUTINE  UPDATES  THE  FILTER  TITB  TRE  MEASUREMENT. 

AUTHOR  :  EARL  MENTZ  AND  PAUL  0.  FILIOS 
VERSION:  2.04 
DATE  :  •  SEP  SS 

PARAMETARS  PASSED: 

Z:  THE  CURRENT  MEASUREMENT 

XH:  THE  CURRENT  FILTER  STATE  ESTIMATE 

HMAT:  THE  FILTER  MEASUREMENT  MATRIX 

CRMAT:  THE  KALMAN  FILTER  GAIN  MATRIX 


PSEUDOCODE: 

COMPUTE  RESIDUAL 

COMPUTE  CHANGE  IN  STATB  ESTIMATE 
ADD  CHANGE  TO  CURRENT  STATE  ESTIMATE 
END 


■« 


Magill,  D.  T.  "Optimal  Adaptive  Estimation  of  Sample 
Stochastic  Processes,"  IEEE  Transact £A  Aut,£ni  a.t  _i.c 
Control .  Vo  1 .  10,  No.  5  :  434-439,  Oct.  1965. 

Hilborn,  C.  6.,  Jr.  and  Demetrios  6.  Lainiotis, 

"Optimal  Estimation  in  the  Presence  of  Unknown 
Parameters,"  IEEE  Transactions  on  Systems  Science  and 
Cybernet ics .  Vol.  5:  38-43,  Jan.  1969. 

Lainiotis,  Demetrios  G.  "Optimal  Adaptive  Estimation: 
Structure  and  Parameter  Adaptation,"  IEEE  Transactions 
on  Automatic  Control .  Vol  16:  160-170,  Apr.  1971. 

_ .  "Partitioning:  A  unifying  Framework  for 

Adaptive  Systems,  I:  Estimation,"  Proceed  ings  of  .the. 
IEEE.  Vol.  64:  1126-1143,  Aug.  1976. 

Maybeck,  Peter  S.  S  tocha  s  tic  Models.  Es  t imat  ion,  and 
Control ,  Vo  lame  2.  New  York:  Academic  Press,  1982. 

Rawkes,  Rechard  M.  and  John  B.  Moore.  "Performance 
Bounds  for  Adaptive  Estimation,"  Proceedings  o f  the 
IEEE.  Vol.  64:  1143-1150,  Aug.  1976. 

Dasgupta,  S.  and  L.  C.  Westphal.  "Convergence  of 
Partitioned  Adaptive  Filters  for  Systems  with  Unknown 
Biases."  IEEE  Transact iops  fin  Automatic  Control . 

Vol.  28:  614-615.  May  1983. 

Athans,  M.,  and  C.  B.  Chang,  "Adaptive  Estimation  and 
Parameter  I d e n t  i  f  i  c a i on  Using  Multiple  Model  Estimation 
Algorithm,"  Technical  Note  1976-28,  ESD-TR-76-1 84 , 

Lincoln  Laboratory,  Lexington,  Mass.,  June  1976. 

Chang,  C.  B.,  and  J.  A.  Tabaczynski,  "Application  of 
State  Estimation  to  Target  Tracking,"  i  EEE 
Transactions  Automat  ic  Control.  Vol.  29,  No.  2:  98-109, 

Feb.  1984. 

Chang,  C.  B.,  and  M.  Athans,  "State  Estimation  for 
Discrete  Systems  with  Switching  Parameters."  IEEE 
Transac  t  ions  on  Ae r oipape  andE.lec_tron.ic  Systems, 

Vol.  14,  No.  3:  418-424,  May  1978. 

Athans,  Michael  et  al.  "The  Stochastic  Control  of  the 
F-8C  Aircraft  Using  a  Mutiple  Model  Adaptive  Control 
(MMAC)  Method  -  Part  1:  Equilibrium  Flight,"  I.EEE 
IXlflAAfiAifiSA  £  A  Auio  m  aii£  Cfint££.l,  Vol.  22,  No.  5:  768- 

780.  Oct.  1977. 

Maybeck,  Pete  S.  EfiiiaAiifiA.  AA4 

£&Al££l>  YfilllBA  i*  ^ew  York:  Academic  Press,  1982. 


Thorp,  James  S.  "Optimal  Tracking  of  Manenvering 
Targets,"  IEEE  Ti42.iicii.0as.  on  A4r0_spa.ee.  and 
Electronic  Si.iie.ma ,  Vol.  9:  512-519,  July  1973. 

Moose,  Richard  L.  "An  Adaptive  State  Estimation 
Solution  to  the  Maneuvering  Target  Problem,"  IEEE 
Transact  ions  fin  Automatic  Control.  Vol.  20:  359-362, 

June  1975. 

Gholson,  N  H.  and  R  L.  Moose.  "Maneuvering  Target 
Tracking  Using  Adaptive  State  Estimation,"  IEEE 
Transac  t ions  on  Aerospace  and  Electronic  Syi.te.m4,  Vol. 
13:  310-317,  May  1977. 

Moose,  Richard  L.  et  al.  "Modeling  and  Estimation  for 
Tracking  Maneuvering  Targets,"  IEEE  Tr.  lane.  Hoai  on 
Aerospace  and  EI44t.ro  nil  Sy4ie.ms.,  Vol.  15:  448-456, 

May  1979. 

Korn,  J.  and  L.  Beean.  Ann  1 i c a t ion  o f  Multiple  Model 
Ada^iiye  Es.i_ima.ii4n  Aijj.or.iihm 4  .to  Maneuver  De.iec.iioa 
and  Estimation.  Contract  DAAK10-82-C-0020.  ALPHATECH, 
Inc.,  Burlington,  Mass.,  June  1983  ( AD-B07 5 92 1 ) . 

Brown,  Grover  R.  "A  New  Look  at  the  Magill  Adaptive 
Filter  as  a  Practical  Means  of  Multiple  Hypothesis 
Testing,"  IEEE  Transact  ions  on  Circuits  and  Systems, 
Vol.  30:  765-768,  Oct.  1983. 

Willsky,  Alan  S.  et  al.  "Dynamic  Model-Based 
Techiques  for  the  Detection  of  Incidents  on  Freeways," 
IEEE  Transact  ions  on  Automat  ic  Coniroi, 

Vol.  24:  347-359,  June  1980. 

Sims,  Craig  S.  and  M.  R.  D'Mello.  "Adaptive 
Deconvolution  of  Seismic  Signals,"  IEEE  1 4441 44  t_i2ni 
on  Geoiiienice  Ele £i r 0 aiaa ,  Vol.  16:  99-103,  Apr.  1978. 

Hostetler,  Larry  D.  and  Ronald  D.  Andreas.  "Nonlinear 
Kalman  Filtering  Techiques  for  Te r r a  in- A ided 
Navigation,"  IEEE  Tr4n4a4t.i2.n4  on  Aut.oma.Lic  Co2_tr.4l, 
Vol.  28:  315-323,  March  1983. 

Mealy,  Gregory  L.  and  Wang  Tang.  "Application  of 
Multiple  Model  Estimation  to  a  Recusive  Terrain  Height 
Correlation  System,"  IEEE  IxiaiJ^iiafil  on  Automatic 
Control  ,  Vol.  28  :  315-323  . 

Moose,  R.  L.,  and  P.  P.  Wang.  "An  Adaptive  Estimator 
with  Learning  for  a  Plant  Containing  Semi-Markov 
Switching  Parameters,"  IEEE  Taaail  c.Li2fil  4S  SyiiES 
Science  a  ad  Cyfc  e  jr  aiiill  •  277-281,  May  1973. 


Meer,  D.  E.  Multiple  Model  Adaptive  Estimation  for 
Space-Time  Point  Process  Observations.  PhD 
dissertation.  School  of  Engineering.  Air  Force 
Institute  of  Technology,  W r i gh t-Pa t t e r s on  AFB,  Ohio, 
Sept.  1982. 

Meer,  D.  E.  and  P.  S.  Maybeck.  "Multiple  Model 
Adaptive  Estimation  for  Space-Time  Point  Process 
Observation,"  Proceedings  pi  p]ip  i EEE  Cppie.xe.iipp  pp 
De c i s ion  api  Con t r o  1  .  Las  Vegas,  Nevada,  Dec.  1984. 

Fry,  C.  M.  and  A.  P.  Sage.  "On  Hierarchical  Structure 
Adaptation  and  Systems  Identification,"  I.p.tprpa.t.ipnp.1 
Journal  pi  Contro 1 ,  Vol.  20:  433-452,  1974. 

Lamb,  P.  R.  and  L.  C.  Westphal.  "  S  i  m  p  1  e  x- d  i  r  e  c  t  e  d 
Partitioned  Adaptive  Filters,"  International  Journal 
of  Control .  Vol.  30:  617-627,  1979. 

Hentz,  K.  P.,  "Feasibility  Analysis  of  Moving  Bank 
Multiple  Model  Adaptive  Estimation  and  Control 
Algorithms."  M.S.  thesis.  School  of  Engineering,  Air 
Force  Institute  of  Technology,  W r  i  ght-Pa 1 1 e r s on  AFB, 
Ohio ,  Dec .  1984 . 

Maybeck,  P.  S.  and  K.  P.  Hentz.  "Investigation  of 
Moving-Bank  Multiple  Model  adaptive  Algorithms," 
Proceedings  pi  the  24 _t£  IEEE  Copierpppp  on  De.pis.ipp  and 
Control.  Ft.  Lauderdale,  Florida,  Dec.  1985. 

Maybeck,  P.  S.  Stochastic  Models.  Est imat  ion,  apd 
Control ,  Vpipme.  3..  New  York:  Academic  Press,  1982. 

Junkins,  J.  L.,  Bodden,  D.  S.,  and  Kamat,  M.  P.  "An 
Eigenvalue  Optimization  Approach  For  Feedback  Control 
of  flexible  Structures."  Proceedings  pi  .the 
Southeastern  Copiprenpp  on  Thepipl  i.pa.1  apd  App lied 
Mecjhanips .  Vol.  II:  303-308,  May  1984. 

Mucken t ha  1 e r ,  T.  V.  "Incorporating  Control  Into  the 
Optimal  Structural  Design  of  Large  Flexible  Space 
Structures"  M.S.  thesis,  School  of  Engineering,  Air 
Force  Institure  of  Technology,  W r  i  gh t-P a 1 1 e r s on  AFB, 
Ohio,  Dec.  1984. 

Jensen,  R.  L.  and  D.  A.  Harnly.  "An  Adaptive 
Distributed-Measurenent  Extended  Kalman  Filter  For  A 
Short  Range  Tracker,"  M.S.  thesis.  School  of 
Engineering,  Air  Force  Institute  of  Technology,  Wright- 
Patterson  AFB,  Ohio,  Dec.  1979. 


Captain  Paul  G.  Filios  was  born  on  June  4,  1957,  in 
Northampton  Massachusetts.  He  graduated  from  high  school  in 
Amherst,  Massachusetts  in  1975.  He  then  matriculated  at  the 
University  of  Massachusetts,  Amherst  campus.  He  graduated 
1979,  with  the  degree  of  Bachelor  of  Science  in  Electrical 
Engineering  and  a  commission  in  the  USAF  through  ROTC  in 
September  1979.  He  was  assigned  to  the  Air  Logistics  Center 
in  San  Antonio  Texas,  Directorate  of  Material  Management; 
and  to  Headquarters,  US  Air  Forces  Europe,  Ramstein  AB  Vest 
Germany,  Deputate  for  Logistics;  before  being  assigned  to 
the  Air  Force  Institute  of  Technology,  in  June  1984. 


Permanent  address:  69  Pelham  Rd 


Amherst,  MA  01002 


r 


I! 

f 


S 


X 

r 


f 

a 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


W  n 


/'♦la.  REPORT  SECURITY  CLASSIFICATION 

'<  UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 

”” lib.  RESTRICTIVE  MARKINGS 


2a  SECURITY  CLASSIFICATION  AUTHORITY 

3b.  OE  CLASSI  FI  CAT  I ON/OOWN  GRADING  SCHEDULE 


3.  OISTRI BUTION/AVAI  LABILITY  OF  REPORT 


A  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

AFIT/GE/ENG/85D-14 


Approved  for  public  release; 
_ distribution  nnl  imi  1-oH 

S.  MONITORING  ORGANIZATION  REPORT  NUMBERIS) 


1 6*.  NAME  OF  PERFORMING  ORGANIZATION  |6b.  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 
I  I  (If  applicable)  1 


School  of  Engineering  |  AFIT/ENG 

6c.  ADDRESS  (City,  State  and  ZIP  Code) 


1 7b.  ADDRESS  (City,  State  and  ZIP  Code) 


Air  Force  Institute  Technology 
Wright-Patterson  AFB,  Ohio  45433 


IBa.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


18b.  OFFICE  SYMBOL  IB.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
inapplicable)  1 


1 8c.  ADDRESS  (City.  Slate  and  ZIP  Code) 


IIP.  SOURCE  OF  FUNDING  NOS. 


PROGRAM 
ELEMENT  NO. 


1 11.  TITLE  (Include  Security  Claetlflcatlon) 


PROJECT 

TASK 

NO. 

NO. 

WORK  UNIT 
NO. 


1 12.  PERSONAL  AUTHOR(S) 

%  Paul  G.  Filios,  B.S.E.E.,  Capt. 

13a.  TYPE  OF  REPORT  13b.  TIME  COVERED 

MS  Thesis  fRQM  to. 

IS.  SUPPLEMENTARY  NOTATION 


USAF 


14.  DATE  OF  REPORT  (Yr..  Mo..  Day) 

1985  Dec 


15.  PAGE  COUNT 

17? 


17.  COSATI  CODES  18.  SUBJECT  TIsRMS  (Continue  on  rfiwnc  if  neeeeeory  and  identify  by  bloc*  number; 

field  group _ sub. gr. _ Adoptive  Control  Systems,  Adaptive  Filters, 

09 _ 05 _ Multiple  Model  Adaptive  Estimation,  Multiple  Model 

lAdaptive  Control _ 

IB.  ABSTRACT  (Continue  on  rwvene  if  neceteary  and  Identify  by  block  number) 


Title;  MOVING-BANK  MULTIPLE  MODEL  ADAPTIVE 
ALGORITHMS  APPLIED  TO  FLEXIBLE 
SPACECRAFT  CONTROL 


Thesis  Chairman;  Peter  S.  Maybeck,  Professor  of  Electrical  Engineering 


KT  C W ^  “w  "  uo-* 

*  /(>  Ohv  Ft 

Iw  Baaaaich  oM  Prnlaaalmiul  Davatapwat 
s*r  Faaea  laatttul*  at  TKhnolopy  (SB6» 
VHpkLMtanaa  AFB  OB  «W3 


1 30.  OISTRI  BUTION/AVAI  LABILITY  OF  ABSTRACT 


21.  ABSTRACT  SECURITY  CLASSIFICATION 


A  m*  _  | 

•  'sunclassifieo/unlimiteo  3  SAMEASRPT.  □  OTIC  USERS  □  UNCLASSIFIED 


1 22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

22b.  TELEPHONE  NUMBER 
(Include  Area  Code) 

32c.  OFFICE  SYMBOL 

[peter  S.  Maybeck,  Prof,  of  Elec.  Eng. 

(513) 255-2057 

AFIT/ENG 

UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGI. 


DO  FORM  1473, 83  APR  edition  of  i  jan  73  is  obsolete.  UNCLASSIFIED _ 

security  classification  of  this  pagi. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  FAGE 


This  investigation  applies  noving-bank  multiple  nodel 
adaptive  e s t iaat ion/ control  algorithas  to  the  control  of  a 
realistic  nodel  of  a  large  flexible  spacecraft.  Moving-bank 
anltiple  nodel  adaptive  estiaation  differs  froa  conventional 
(fnll-bank)  anltiple  aodel  adaptive  satiation  in  that  a 
substantially  reduced  nnaber  of  elemental  filters  is 
required  for  the  noving-bank  eatiaator  (9  va.  100  for  the 
system  modeled  in  this  thesis).  The  positions  in  paraaeter 
space  that  the  reduced  nuaber  of  eleaental  filters  occupy 
are  dynamically  re-declared;  i.e..  the  aoving-bank  slides 
about  the  parameter  space  in  search  of  the  true  parameter 
point. 

Critical  to  the  performance  of  the  noving-bank  anltiple 
aodel  adaptive  estimator  is  the  decision  logic  used  to 
determine  which  eleaental  filters  are  iapleaented  in  the 
bank,  and  when  to  change  this  decision.  The  decision  logics 
discussed  focus  on  three  situations:  initial  acquisition  of 
the  unknown  paraaeter  values,  through  reducing  bank 
discretization;  tracking  the  unknown  paraaeter  valnes, 
through  bank  aovement;  and  reaoquiaition  of  the  unknown 
parameters  following  a  large  juap  change  in  their  values, 
through  expanding  bank  discretization.  Aabiguity  function 
analysis  is  used  to  predict  performance  in  these  situations. 

The  systea  to  be  controlled  is  a  siaplif  led  model  of  a 
large  scale  space  structure.  Its  equations  of  notion  are 
developed  and  placed  in  state  space  fora,  the  states  being 
the  positions  and  velocities  of  the  rigid  body  mode  and  the 
second  and  fourth  bending  nodes.  The  state  space  matrices 
describing  the  systea  are  computed  based  on  noninal  values 
for  all  physical  parameters  with  the  exception  of  the  aass 
density  of  the  structure  arns  and  their  aodulus  of 
elasticity.  These  two  paraaeters  are  allowed  to  vary  in 
discrete  steps,  estaplishing  the  paraaeter  space.  It  is 
then  attempted  to  control  the  states  to  the  quiesent  state, 
using  moving-bank  multiple  aodel  adaptive  algorithas. 

The  results  indicate  that,  although  the  systea  under 
study  did  not  have  a  great  need  for  adaptive  estiaation  and 
control,  the  anltiple  aodel  adaptive  estiaator  performs 
essentially  identically  to  a  single  filter  artificially 
knowledgeable  of  the  uncertain  paraaeter  values.  In 
addition  changing  bank  discretization  for  the  initial 
parameter  aquisition  phase  speeded  acquisition.  However, 
the  bank  was  unable  to  expand  following  a  juap  change  in  the 
uncertain  paraaeter  values,  in  order  to  restart  the 
acquisition  phasep  the  bank  tracked  the  juap  change  through 
movement  alone.  Aabiguity  function  analysis  proved  to  be  an 
excellent  predictor  of  bank  performance,  and  should  be  used 
as  a  design  tool. 
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